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ABSTRACT 

We compare the yields of "^"^Ti and ^^Ni produced from post-processing the ther- 
modynamic trajectories from three different core-collapse models - a Cassiopeia A 
progenitor, a double shock hypernova progenitor, and a rotating 2D explosion - with 
the yields from exponential and power-law trajectories. The peak temperatures and 
densities achieved in these core-collapse models span several of the distinct nucle- 
osynthesis regions we identify, resulting in different trends in the ''^Ti and ^^Ni yields 
for different mass elements. The "^Ti and ^^Ni mass fraction profiles from the expo- 
nential and power-law profiles generally explain the tendencies of the post-processed 
yields, depending on which regions are traversed by the model. We find integrated 
yields of "^"^Ti and ^^Ni from the exponential and power-law trajectories are generally 
within a factor 2 or less of the post-process yields. We also analyze the influence of 
specific nuclear reactions on the '^'^Ti and ^^Ni abundance evolution. Reactions that 
affect all yields globally are the 3a, p(e~, Ve)n and n(e^, v^)p. The rest of the reactions 
are ranked according to their degree of impact on the synthesis of '*'*Ti. The primary 
ones include ^Ti(a, pY^Y, ^^Ca(a, y)^Ti, ^^Y(p, yY^Cr, "^^Caia, pY^Sc, ^^F(a, pf^Ne, 
^^Na(a,p)^Mg, ^^Sc(p,yy^Ti, ^^Sc(p,y)^Ti, ^^Ti(p,yy^Y, and 5^Ni(p, 7)5^Cu, along 
with numerous weak reactions. Our analysis suggests that not all '^'^Ti need be pro- 
duced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria 
in combination with timescale effects for the expansion profile may account for the 
paucity of "^"^Ti observed in supernovae remnants. 
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1. Introduction 



Core-collapse supernovae inject energy and material enriched freshly synthesized isotopes 
into the interstellar medium. Some of this material is short to medium lived radioactivities with 
half-lives ranging from several days to several million years. Detecting y-rays from the decay 
chains of such isotopes, either in individual supernova remnants or through the accumulation 
of material in interstellar medium, provides a direct calibration of the nucleosynthesis in core- 
collapse events. For example, the radioactive decay '^'^Ti and ^^Ni has significant observational 



consequences for the light curves of core-collapse supernovae ( [Arnett et al.||1989[ [Timmes et al. 
T996| IVink et al.||200T| [Hungerford eTaL][2005| [Renaud et al.][2006t [Young et al.|[2006l ), isotopic 
patterns measured in primitive meteorites ( Wadhwa et al.|2007 1 and presolar grains ( Zinner|1998 ), 
anomalies in a deep-sea crust ( Knie et aLl|2004 ), and the solar abundances of ^"^Ca and ^^Fe (Lod- 
ders|2003D . 



The past decade has brought substantial progress to the theory of core-collapse supernovae. 
There now seems to be general agreement that hydrodynamic instabilities above the proto-neutron 
star play a crucial role in, not only achieving an explosion, but also in determining critical prop- 



erties such as the timing, strength and asymmetry of this explosion ( [Buras et aL||2006[ |Bruenn 
eraL][2006| [Kifonidis et aL][2006t [Fryer & Youngl[2007| [Messer et al.][2008t [Ott et al.[[2008| |Lu 
nardini et"aL][2008[ ). For explosion scenarios where the growth of these instabilities is sufficiently 
long, neutrino transport through this region seems capable of resetting the electron fraction of 
at least some material from being neutron-rich to being proton-rich ( Pruet et al.|2005[ 2006 , Buras 
et al.[[200^ ; Frohlich et al.[[2006 ). While the details are sensitive to the numerical techniques and 
physical approximations employed in different simulations, the range of explosion strengths and 
timings obtained imply significant variations in the evolutions of the temperature, density and 
in the tumultuous inner regions. 

Observations of '^'^Ti and ^^Ni in individual core-collapse supernova may provide the best 
probes for constraining aspects of the explosion mechanism precisely because the production of 
these two isotopes are sensitive to the temperature, density and Y^ evolution. Perhaps most com- 
pelling are abundance determinations of the Cassiopeia A remnant from Compton Gamma Ray 
Observatory, BeppoSAX, INTEGRAL, and Chandra measurements. The inferred ratio of "^"^Ti to 
^^Ni in Cas A is higher than that predicted by standard, spherical supernova explosion models 



(Young et al. 2006, Young & Fryer 2007). The solar abundance ratio of '^'^Ca to ^^Fe is similar 



to Cas A's ratio of "^^^Ti to ^^Ni, suggesting that spherical models are simply falling short in their 
synthesis of "^"^Ti. Of course, it can be argued that Cas A was simply a peculiar event (The et al 



20061 ). 



Multi-dimensional effects may play some role in resolving this discrepancy (e.g., Arnett et al. 
2008[ ). Explosions with artificially imparted asymmetries in 2D were modeled by Nagataki et al. 
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( 1997 ) to show that bipolar explosion scenarios could account for enhanced "^"^Ti synthesis along the 
poles of model supernova explosions. Simulations of core-collapse and hypernovae, where high 
energies and large asymmetries are imparted to launch the explosion, can reproduce the trends in 



the abundances of metal poor stars and imply larger masses of Ti are ejected (Tominaga et al. 



2007, Umeda & Nomoto 2008). If effects from asymmetries are important for setting the nu- 
cleosynthesis of '^Ti and ^^Ni, then quantifying the physical and numerical uncertainties which 
determine those asymmetries becomes important. 

However, multi-dimensional explosion simulations are resource intensive, and thus run pri- 
marily to address hydrodynamic and transport aspects and uncertainties of supernovae. Such mod- 
els have not yet been run to assess the sensitivity of isotopic yields ( [Young & Fryer|[2"007 l) to the 
nuclear physics input. Parameterized expansion profiles bypass these difficulties by simplifying 
the hydrodynamics in favor of focusing on nucleosynthesis. A motivation for this paper is to begin 
the process of examining the interplay between these two modes of analysis. Thus, in this paper we 
focus on the production of '^'^Ti and ^^Ni from classic adiabatic freeze-out thermodynamic trajecto- 
ries, power-law thermodynamic trajectories suggested by 2D explosion models, and core-collapse 
supernova models. We explore in detail the sensitivity of the "^^^Ti and ^^Ni produced to variations in 
the reaction rates, electron fraction, and nuclear network size with the simple thermodynamic tra- 
jectories. We assess how yields determined from the simple thermodynamic trajectories compare 
to the post-process yields from complex simulations of core-collapse supernovae. This assessment 
offers a calibration of where simple trajectories provide a reasonable approximation to the final 
yields, and allows discovery of which regions in the explosion models deviate from the simple 
trajectories and why they differ. Previous efforts along these lines explored the sensitivity of ^"^Ti 



synthesis to the assumed reaction rates or the electron fraction (Woosley et al. 1973, Woosley & 



Hofi^man|199"2HThe et al.|1998[ [Hoffman et al.|2010| ). In this paper we study the sensitivity of ^"^Ti 



and ^^Ni synthesis for both dependencies over an extended parameter space. 

In §|2]we briefly discuss equilibrium states and in ^we present the exponential and power-law 
thermodynamic trajectories to be interrogated. Section |4] considers general trends of "^"^Ti and ^^Ni 
from these trajectories in the peak temperature-density plane. We also show where in this plane 
multi-dimensional models of asymmetric supernovae and hypernovae tend to reside. In ^ we 
discuss the nucleosynthesis of '^'^Ti and ^^Ni in material with different Y^, while the sensitivities to 
reaction rate values and network size are discussed in §|6]and ^respectively. Section [8] describes 
the yields of "^"^Ti and ^^Ni from post-processing core-collapse trajectories, compared with the 
yields from the parameterized profiles. We conclude with a summary of our main results in §|9j 

We establish our nomenclature and conventions. Let isotope i have Z, protons. A, nucleons 
(protons -I- neutrons), and an atomic weight Wj. We shall assume Wj = Ai. Let the aggregate 
total of isotope / have a baryon number density (in cm'^) in material with a temperature T (in 
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K) and a baryon mass density p (in g cm"^). Define the dimensionless mass fraction of isotope 
i as Xi = pilp = HiAilipNA), where Na is the Avogadro's number, and the molar fraction of 
isotope / as 7, = Xi/Ai. The electron fraction, or more properly, the total proton to nucleon ratio is 
Ye = Yji^iYi = Tji^i/^i ^i- Wc define "nuclear flow" to mean the instantaneous rate of change of 
isotope i's molar abundance with time, dYi/dt, due to a given nuclear reaction ( Iliadis||2007| ). For 
any single reaction linking isotope / with isotope j there is a forward flow, a reverse flow, and a 
relative net flow 0,=(forward - reverse)/max(forward,reverse) that measures the equilibrium state 
of the reaction. 



2. Silicon burning and equilibrium states 

Silicon burning is the last exothermic burning stage and produces the Fe-peak nuclei. Due 
to Coulomb repulsion, it is rather improbable that two ^^Si nuclei will fuse to ^^Ni. Instead, a 
photodisintegration driven rearrangement of the abundances takes place, originating from equilib- 
ria established among individual reactions with their reverses ( [Bodansky et al.||1968 1. When such 



equilibria happen among many reactions, the plasma reaches an equilibrium state where nuclei 
merge into clusters. Units of interaction are no longer nuclei, but the clusters themselves, which 
adapt their properties according to the local thermodynamic conditions. In general, not all reac- 
tions are in equilibrium. Consequently, this state is named quasi-static equilibrium (henceforth 
QSE). The special case where all strong and electromagnetic reactions are balanced by their re- 
verses is called nuclear statistical equilibrium (henceforth NSE), because all mass fractions may 
be described in terms of statistical properties of excited nuclear states (partition functions) and nu- 
clear structure variables (masses and Q values). Weak interactions are always excluded from these 
definitions, since for conditions relevant to hadronic physics they never attain equilibrium. Hence, 
equilibrium notions are related only with strong and electromagnetic interactions. In practice,there 
is either one cluster in NSE or QSE, or two QSE clusters, one for the Si-group and one for the 
Fe-group nuclei. 

An NSE state may be completely described by a triplet of macroscopic parameters such as 
temperature, density and electron fraction Y^. A QSE state requires additional parameters, one 
for each equilibrium cluster, which may be chosen to be the number of nuclei in each cluster 
( [Meyer et al.||1998| Wallerstein et al.||1997 ). The mass fractions of nuclei in such equilibria states 



are completely described as functions of these parameters. A benefit from this property, is that the 
choice of initial composition has no impact on the equilibrium state, as long as it remains consistent 
with the equilibria parameters. Thus, if an equilibrium state is established, the details how the 
plasma attained that equilibrium are not necessary to model aspects of the continuing evolution. 
This feature is the basis for reliable results from parameterized expansion profiles, whose starting 
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point is the moment the explosion shock strikes the inner stellar layers. 

Reaction cross sections have, in general, an asymptotic trend towards a saturation value at 
high energies. Reactions tend to get balanced by their reverses in this regime. A reaction rate is 



connected to its reverse according to the detailed balance theorem ( |Iliadis||2007] ). The dominant 
term in this relationship is exp{-Q/kT), implying that the reaction Q value and the temperature are 
the foremost magnitudes related with the ability of a reaction to reach and maintain equilibrium. 
There is a linear density dependence to this relationship only for reactions involving photons. 
Large Q values result in sensitive reaction equilibria, which are the first to break for decreasing 
temperature. However, these reactions become the most efficient flow carriers once they break 
equilibrium, since they release the largest amounts of energy per reaction. Figure [T] shows the 
Q value distribution for alpha particle captures within our base network containing 489 isotopes, 
which includes all reactions that may directly affect '^'^Ti. Table [Tjgives a complete specification of 
all networks used in this study. 

The temperature expresses the internal energies of the nuclei, while the density is related to 
their availability for reactions. As a result, the larger the temperatures and densities are, the more 
equilibria exist. This effect allows NSE to be established for large temperatures and densities. 
For smaller temperatures and densities, certain equilibria start breaking, but a large scale QSE 
structure in the plasma may still exist. For yet lower thermodynamic conditions, the large scale 
clusters dissolve into smaller clusters. Very low temperatures and densities are not adequate to 
establish a significant amount of equilibria. A few isolated equilibria may exist, but without any 
specific connection between them. Depending on the initial peak temperature and density, the 
plasma may experience one or more of these states during the expansion of the ejecta (see Figure 
[2]). Unfortunately, the threshold conditions to border each regime cannot be known accurately, 
since they are sensitive to the number of species involved (network size), the Q values of the 
associated reactions and the reaction rates used. Such borders exist in nature though, and the time 
spent by the plasma in each regime may affect the final yields. 

External flow supply to reactions may also result in equilibria breaks, even for constant tem- 
perature and density. When there are no external abundance flows to reaction k)m, the condition 
0( = -0m nieans the forward flow for isotope / is the reverse flow for isotope m and vice versa. 
Due to this condition, this reaction may achieve equilibrium, a case where both and (pm are equal 
to zero. Assume now an external abundance flow, say from another reaction that is not in equilib- 
rium, that supplies flow only to isotope /. As long as the external flow is significant in magnitude, 
^ I0ml because the external flow is additive only to 0,-. That is, the external flow term is applied 
to the equation for dYt/dt but not to the equation for dY,„/dt, and the reaction z(j, k)m breaks equi- 
librium. If the external abundance flow is applied long enough, a new equilibrium state may be 
established. Turning off the external flow causes the reaction to be driven back to equilibrium with 



-6- 



a new abundance ratio, depending on the flow transfer by the external agent. In general, starting 
from an equilibrium state, a transition to any other state means some \(()\'s must be greater than other 
101 's during the transition. External flows may also affect the equilibrium state of small clusters of 
nuclei, or even large scale QSE clusters (see 



3. Parameterized Thermodynamic Trajectories 

We use two parameterized expansion profiles to identify robust trends and uncertainties in the 
yield of "^"^Ti and ^^Ni. Both profiles assume that a passing shock wave heats material to a peak 
temperature Tq and compresses the material to a peak density po- This material then expands and 
cools down (freezes out) under the assumption of a constant T^/p evolution (radiation entropy in 
suitable limits) until the temperature and density are reduced to the extent that nuclear reactions 
cease. Our adiabatic freeze-out trajectories ( |Hoyle et al.|1964 ; Fowler & Hoyle|1964 ) 



dt 3t dt T 
T(t) = Toexp(-t/3T) p(t)=poexp(-t/T) (2) 
are used with a static free-fall timescale for the expanding ejecta 

T = (lAnGpo)-"^ « 446/p;/' s (3) 

In this formulation the temperature and density evolutions are decoupled. If one uses p(t) in the 
expansion timescale instead of the peak density po, then the temperature ordinary differential equa- 
tion becomes coupled to the density evolution. 

The second thermodynamic profile we use is based on homologous expansion. For a fixed 
expansion velocity v, the distances increase as R{t) = Ro + vt, the density scales as p{t) ~ 1 /Rit)^ = 
1 /{Rq + vtf and the temperature scales through p cc T^. Specifically we use 

dT _ -To dp _ -3po 

dt ~ 1/2(2?+ 1)2 dt ~ l/2(2f+ 1)4 

r(,) = ^ P(0 = (5^, (5) 

where the coefficients in the denominator are chosen to mimic trajectories taken from core-collapse 
simulations. Substituting the power-law solution into the ordinary differential equations they orig- 
inate from and eliminating the direct time dependence 

dT _ Tq dp _ 3po 

dt ~ 0.5(p/po)2/3 dt ~ 0.5{T/Toy 
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shows the temperature and density evolutions are coupled for the power- law trajectories. 

Figure |2] compares the general properties of these two parameterized profiles. For a given 
initial condition, the power-law evolution is always slower than the exponential one. Moreover, 
the power-law evolution becomes slower for increasing initial values. The differences in these two 
profiles affect the final yields of "^"^Ti and ^^Ni as material traverses different burning regimes on 
different timescales. The figure also depicts the NSE, global QSE, local QSE, and final freeze- 
out burning regimes. The exponential and power-law trajectories are chosen so that they bound 
in general the temperature and density trajectories of particles from the '^'^Ti and ^^Ni producing 
regions of spherically symmetric and 2D explosion models. 

For any given peak temperature, peak density, and initial electron fraction Yg we want to know 
the mass fraction of "^'^Ti and ^^Ni produced by nuclear burning from the exponential and power- 
law profiles. We chose peak temperatures, peak densities, and values spanning the range of 
4 X 10*^ < To < 10 X 10'' K, 10^ <po< 10^° g cm-\ and 0.48 < 7, < 0.52. This parameter space 
covers the conditions encountered in most core-collapse supernova models which produce some 
'^'^Ti or ^^Ni. When sampling this parameter space between these limits we use 121 points, equally 
spaced in base 10 logarithm, for the peak temperature or density and increments of 0.002 in Yg. 
That is, for any value of Y^, we compute the final nucleosynthesis at 121x121 points in the peak 
temperature-density plane using mature reaction networks ( [Timmes 1999[ Fryxell et al. 2000). 



Using a larger number of sample points does not alter our main results and conclusions. Our initial 
composition for any starting (Tq, po, Y^) triplet is pure ^*^Si for symmetric matter {Ye=0.5). We then 
added protons or neutrons to make initial composition either proton or neutron-rich respectively. 
Specifically, we used X(^*^Si) = 1 - 12 7^. - 1| and either X(p) = |2 7^. - 1| for proton-rich compositions 
(Ye > 0.5) or X(n) = |2 7^ - 1| for neutron-rich compositions (7g < 0.5) to set the initial ^^Si, proton 
or neutron mass fractions. As we show below in ^ the choice of ^^Si is not important for vast 
regions of the chosen parameter space. 



4. Trends in the Peak Temperature-Density Plane 

Figure [b] shows the mass fractions of '^'^Ti, ^^Ni and '^He produced in the peak temperature- 
density plane for the exponential and power-law profiles and an initial 7^=0.5. Each point in 
the plane represents one set of initial conditions, which are evolved forward in time according to 
Equations [T] or |4| with the final freeze-out abundance of "^"^Ti and ^^Ni recorded. The color map 
is logarithmic, spanning mass fractions from 10"^ to 10 for "^"^Ti and from 1 to 10"^° for ^^Ni. 
The overlaid colored triangles correspond to the temperature and density of particles from a suite 
of supernovae and coUapsar simulations in the region where '^'^Ti and ^^Ni are produced. Not all 
particles have an initial 7^=0.5, but are relatively close to it. Each supernova model generally spans 
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the full range of peak temperature, but the peak density is confined to a strip of one or two orders 
of magnitude. 

Several striking patterns emerge from these contour plots. The first is "^"^Ti seems to be pro- 
duced overall with an average mass fraction X(^Ti) ~ 10"^, except in certain regions where it 
gets depleted. The depletion region extends along a thin line for low temperatures (oriented ap- 
proximately 70° with respect to the temperature axis), smoothly bending over into a wider, more 
horizontal band for relatively high temperatures and densities. We name this depletion region the 
"chasm". The chasm separates the peak temperature-density plane into distinct regions controlled 
by different burning processes. 

The second pattern is that the "^"^Ti contour plots for the two thermodynamic profiles have 
the same general structure, except that the chasm for the power-law profile is located at lower 
densities and slightly wider compared to the exponential profile. Hence, the power-law chasm 
begins to encompass the majority of the overlaid particles. It is possible that the existence of the 
chasm region is why so few supernova remnants have been observed in the glow of radioactive 
"^"^Ti. The total mass of "^Ti ejected by an individual core-collapse supernova depends critically on 
(i) the location of its thermodynamic points in the peak temperature-density plane and (ii) the exact 
expansion profile that the ejecta follow past the explosion. The impact of the latter is expressed 
as the chasm's ability to "shift" and "widen" itself from exact profile to exact profile. The third 
pattern is ^^Ni has large mass fractions and is relatively featureless in the peak temperature-density 
plane. Large variations in observed '^'^Ti to ^^Ni ratios are primarily due to variations in '^'^Ti. 

The chasm's formation and trends with thermodynamic history are the primary motivation 
for using two parameterized profiles. Our analysis to uncover the nuclear physics controlling the 
chasm is two-fold. First, we ascertain the basic synthesis mechanisms of '^'^Ti in distinctive ther- 
modynamic regions through a series of nuclear reaction network calculations (section [5]). Second, 
we identify reactions crucial to '^"^Ti in each thermodynamic region via a three-stage process based 
roughly on the methodology established by The et al. (1998]), but modified because we are inter- 



ested in more than one normalization point in each peak temperature-density plane (section [6]). 



5. Nucleosynthesis of Ti and ^^Ni 

Figure [4] shows the mass fraction of "^"^Ti in the peak temperature-density plane for the expo- 
nential profile and 7^=0.5. Each point in this plane represents one set of initial conditions which 
evolve forward in time with the final freeze-out abundance of "^'^Ti cataloged. Six different regions 
are labeled which are characterized by specific nuclear burning patterns controlling the produc- 
tion of '^'^Ti. Despite the differences in timescale, the corresponding temperature-density contour 
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plot following the power-law profile contains the same number of regions, with similar physics 
associated between regions of the same label. Thus, the duration of the hydrodynamic expansion 
does not explain the underlying structure of the contour plots for the final ''^'^Ti and ^^Ni yields. 
Instead, the entropy during the expansion drives the nucleosynthesis, by affecting the strengths 
of key nuclear reactions, and causing phase transitions in the burning process for certain critical 
temperatures and densities. The phase transitions are followed by a change in the burning state. 
On the other hand, the expansion timescale affects the locus of the borders among different re- 
gions on the contour plot (Figure |3]). For increasing expansion timescale the plasma spends more 
time within each burning state. Depending on the region in the peak temperature-density plane, 
the evolution may include some or all states between NSE and non-equilibrium nuclear burning 
(see Figure [2]). Timescale differences between profiles result in different density values, when the 
temperature acquires a threshold value indicative of a phase transition. Since both profiles attain 
constant radiation entropy, different densities at threshold temperatures translate to different peak 
densities and hence, border shifts between regions on the temperature-density plane (e.g. see § |6.3| 
for the explanation of the chasm shift). 

For most of the regions in Figure |4| the peak conditions are sufficiently large that the plasma 
reaches a large scale equilibrium state (NSE or QSE) on timescales much shorter than the freeze- 
out timescale. During the first time steps of a reaction network calculation the initial composition 
rearranges to an NSE or a QSE distribution well before the temperature and density begin evolving. 
This rapid rearrangement appears as vertical line in many of our plots. As the plasma subsequently 
cools and rarefies the first transition occurs when the NSE state can no longer be maintained. The 
threshold temperature for NSE is usually taken to be = T K) ~ 5. The density at this 
threshold temperature determines the subsequent burning process by prescribing both the available 
amount of nuclear fuel and the dominant flows that consume the fuel. 

Region 1 is essentially a freeze-out from NSE, henceforth termed a "normal freeze-out" 
dWoosley et al.|1973|peyer|1994|peyer et al.|1998HHix & Thielemann|19"99l ). When the temper- 



ature falls to the = 5 threshold temperature, the density is ~ 1.0x10^ g cm"^ for the high peak 
density region above the horizontal band of the chasm. At this density an NSE distribution is domi- 
nated by ^^Ni, contains a significant amount of Si-group and Fe-group nuclei, but a relatively small 
amount of free alpha particles (Xq, ~ 10"^). This density is large enough to favor particle captures, 
but the temperature is such that photodisintegration reactions are not negligible either. The large 
scale equilibrium structure is maintained until complete freeze-out, since Xq, «; 1 for the majority 
of equilibrium states and the 3a reaction is always dominated by its inverse photodisintegration. 
Since normal freeze-out is a dynamic process though, some individual equilibria are broken as the 
plasma cools and rarefies and QSE estimates become progressively more accurate compared to 
NSE estimates ( Woosley et al.|1973] ). Yields for the isotopes plotted in Figure |5]for region 1 (first 



row of plots) are not far from NSE or QSE yields. Thus, network calculations in this region may 
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be avoided and accurate estimates for yields may be determined only by nuclear properties such as 
masses and Q values. 

Equilibrium estimates for Si-group and Fe-group nuclei during the initial stages of the expan- 
sion remain relevant for region 2. When the temperature falls to the = 5 threshold temperature 
the density is ~ 1.0x10^ g cm"^. The low availability of alpha particles at this density does not al- 
low the 3a reaction to dominate its inverse, preventing significant flow from the light nuclei to the 
equilibrium cluster. Compared to region 1 though, not all of the capture reactions have the same 
efficiency. The large a thresholds in nuclei between N,Z=20 and N,Z=26 closed shell configura- 
tions results in a phase transition which is responsible for the formation of the chasm. Because of 
the large Q values associated with a capture in the mass range 42 < A < 53 due to shell structure 
(Figure [T]), these reactions are the first to break the local equilibria and form a continuous passage 
of nuclear flow from the Si-Ca-group to the Fe-group nuclei. The large equilibrium cluster dis- 
solves into two smaller ones, with ^"^Ti being located within the upper mass limits of the Si-Ca 
cluster, while ^^Ni is centralized in the Fe-group. The flow transfer between the two equilibrium 
clusters results in the depletion of '^'^Ti and the rest of the isotopes in the Si-group by the end of the 
thermodynamic evolution (second row of Figure [s]). On the contrary, ^^Ni is one of the Fe-group 
isotopes that benefit from this transfer since the reaction equilibria in its neighborhood are main- 
tained until freeze-out. Equilibria estimates for this small group of nuclei within the Fe-group are 
still a good approximation. The formation of the chasm in Figure |4] is a direct result of a phase 
transition from the single cluster QSE configuration to a double cluster QSE configuration and the 
subsequent flow leakage. 



Region 3 corresponds to the conditions of a-rich freeze-out (Woosley et al. |1973| ). As the 



plasma cools and rarefies, most Si-group and Fe-group mass fractions acquire the topology of 
an "arc" in going from low values at high temperatures to a local maximum and back to a local 
minimum at cooler temperatures, while in QSE (third row of Figure [5]). The density at the Tg=5 
threshold temperature within region 3 spans 10"^ < p < 10^ g cm^^, resulting in less efficient 
particle captures compared to regions 1 and 2, and a helium mass fraction X„ ~ 10"^ The excess 
of free alpha particles allows the 3a rate to dominate its inverse photodisintegration, leading to a 
new phase transition. Although the 3a rate remains relatively slow ( |The et aL] 1998 1, it supplies 



external flow which breaks the local equilibria in the neighborhood of '^C, ^^O and ^"^Mg. The 
subsequent energy release from alpha capture reactions provides a significant nuclear flow towards 
heavier nuclei by breaking successively other local equilibria. The QSE cluster changes its shape 
and shifts gradually upwards in mass, instead of dissolving into two clusters. Meyer et"aL] ^99S ) 



identified this cluster motion based on QSE calculations. The mass fractions of nuclei which 
suddenly find themselves outside the QSE cluster begin an ascending track. These are primarily 
the Si-group nuclei (including '^'*Ti), and a few from the Fe-group. Close to complete freeze-out, 
the yields for these nuclei are orders of magnitude larger than their corresponding minimum value 
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reached prior to the phase transition. Because the 3a reaction itself is not very efficient, the process 
ends up with an excess of alpha particles Xq, ~ 10"^. 

Large scale QSE calculations cannot account for the increase of the "^Ti mass fraction curve 
past the arc, since "^"^Ti and other related isotopes have decoupled from the large scale equilibrium 
cluster. However, the phase transition is not abrupt in shifting from total equilibrium to total non- 
equilibrium. Nuclear flow analysis shows that ^"^Ti participates in a smaller, local equilibrium 
pattern which is responsible for its ascending trend in the mass fraction curve at the end of freeze- 
out. This transition is demonstrated in Figure [6] which displays the reaction links between f7/2-shell 
nuclei located between the Z,N=20 and Z,N=28 closed shells. The top panel in Figure [6] shows 
the network of reactions prior to the "^"^Ti abundance minimum, which is characterized by (p, y) 
equilibria along the N=22 and N=24 isotone chains connected by (a, p) and (or, y) channels in 
equilibrium. These equilibria guarantee linkage of '^'^Ti to the large scale QSE cluster, and hence, 
the downward portion of the mass fraction curve is produced. The breakdown of the equilibrium 
conditions for the '*'^Ti(Qr, p)'^^V link signals the phase transition for "^"^Ti. Soon afterward, the rest 
of the (a, p) and (a, y) equilibria connecting the N=22 and N=24 isotone chains break, as reflected 
in the increase of actual net flow shown in the middle panel in Figure [6| "^"^Ti is left in (p, y) 
equilibria along the N=22 isotone chain with '^^V, '^^Cr, ^'^Mn and "^^Fe and its mass fraction starts 
to increase from the local minimum (bottom panel of Figure [6]). It is this (p,y) equilibria chain 
which is responsible for the rising portion of the mass fraction curve after the local minimum. This 
pattern persists until complete freeze-out. 

These reaction network flow study results can be verified by localized QSE calculations. The 
advantage of QSE modeling is that the abundances of all isotopes within a cluster may be expressed 
by a semi-analytical formalism, in terms of the network abundances of free protons, neutrons and 



an arbitrarily chosen reference isotope. For this purpose, we adopt the Hix & Thielemann ( 1996 



1999[ ) formalism. We model the cases of equilibrium (i) between the N=22 and N=24 isotone 



chains and (ii) only along the N=22 isotone chain throughout the evolution, corresponding to the 
top and bottom panels in Figure |6] respectively. Both cases reproduce the arc topology of the "^"^Ti 
mass fraction curve. However, the first case does not reproduce the ascending part of the curve 
beyond the local minimum in Figure [v] Instead, the '^'^Ti curve continues to descend, expressing 
the trend of '^'^Ti, were it to remain in global QSE. The second case on the other hand, which 
expresses only (p, y) equilibria, fits the network results until the point where complete freeze- 
out occurs. The discrepancy beyond this point relies on the fact that nuclear reactions no longer 
take place. Thus, mass fractions do not change any more and the curve from network calculations 
acquires a plateau. This general behavior applies to most of the elements within the silicon and iron 
groups, as demonstrated for a small subset within these groups at the upper right panel in Figure|7] 
The mass fraction trends of an isotope depend strongly on the local reaction equilibria within its 
neighborhood. Further equihbria isotone chains are readily identifiable in Figure [6} For example. 
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crucial equilibria reactions for "^"Ca are '*''Ca(a, y)'^Ti and '*^Sc(p, y^'^Ti. Its mass fraction profile 
in the lower left of Figure |7] is in accordance with the general mechanism. The increase of the 
"^^Ca mass fraction with cooling is maintained through the (p, y) equilibria along the N=20 isotone 
chain. A similar development can be observed for "^^Cr as shown in the lower right of Figure [v} the 
crucial reaction now is '^^Cr(a, p)^^Mn. 

Region 4 is a special case of an a-rich freeze-out. Within this region, the p(e~,Ve)n and 
n(e'^, v^)p reactions exert a greater influence compared to the other regions. These reactions drive 
the composition slightly proton-rich near the beginning of the evolution when temperature and 
density are still large. Impacts to the burning processes for proton-rich composition are described 
in more detail in ^5.1[ but some of the impacts include a relatively high number of free protons 
and an enhanced efficiency of proton captures ( |Pruet et al.|2005[ |2006[ |Buras et al.|2006^ [Frohlich 



et al.||2006| ). Consequently, this region is a proton-rich, a-rich freeze-out, henceforth an "or/?- 
rich freeze-out". The fourth row of Figure [s] shows the mass fraction profiles for "^^^Ti in region 
4 have certain similarities to the profiles in region 3. A characteristic arc of large scale QSE is 
formed, followed by the ascending track due to the equilibrium chain connecting "^"^Ti, "^^V, "^^Cr, 
"^^Mn and "^^Fe via (p, y) reactions along the N=22 isotone chain. Among these linking reactions 
'^^V(p, yY^Cr has the largest Q value, and thus will break its equilibrium first as the plasma cools 
and rarefies. When this reaction breaks equilibrium, the N=22 isotone chain dissolves into two 
smaller clusters, the first between "^"^Ti and "^^V and a second between "^^Cr, "^^Mn and "^^Fe. This is 
the second phase transition that '^'^Ti sustains during its evolution. Similar transitions occur along 
other isotone chains. Flows are now carried among such isolated small scale clusters by out of 
equilibrium alpha and proton captures. These flows favor mostly the proton-rich nuclei, resulting 
in a decrease for "^"^Ti and other symmetric isotopes. Thus, a second arc is clearly identifiable in 
the mass fraction curve for most of the isotopes in the fourth row of Figure |5j The ascending 
track beyond the second arc for '^'^Ti and most of the symmetric isotopes is a consequence of the 
flow transfer through weak interactions at the expense of proton-rich nuclei, when the strong and 
electromagnetic reactions become ineffective as freeze-out takes place. 

In region 5 the temperatures are initially large enough to establish equilibrium (NSE or QSE), 
but the initial densities are so low that photodisintegrations soon dominate capture reactions. Long 
before the complete freeze-out, all nuclei dissolve into neutrons, protons and a-particles. A slight 
recombination takes place during the final stages of the freeze-out producing traces of '^C, ^^O 
and ^^Si. The recombination is driven mostly by the 3a, p(e , Ve)n and n(e^,v^)p reactions in 
Yg > 0.5 matter. The products of this recombination set the seed for a following chain of (p, y) 
and (n,p) reactions that produce heavier elements, including "^"^Ti, ^^Ni and the heaviest isotopes 
in the network used for the calculation. Similarly to region 4, weak interactions at the close of the 
process carry some flow from proton-rich nuclei to symmetric ones, enhancing this way the mass 
fractions of '^'^Ti and ^^Ni. However, the contributions of the recombination and the chain of (p,7) 
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and (n, p) reactions are not adequate to yield large production factors for most of the isotopes. The 
final composition is dominated by free alpha particles and protons, establishing this region to be a 
photodisintegration driven regime. 

Region 6 represents incomplete silicon burning, where ^^Si gradually dominates ^^Ni from 
region 1 to the left of the thin chasm line towards the inner part of this region. The peak tempera- 
tures and densities are such that the timescale to reach a single cluster QSE state is comparable or 
larger than the expansion timescale. Multiple small scale QSE clusters are formed, but they do not 
merge successfully into one large scale cluster. The mass fractions freeze out from the established 
equilibrium state without sustaining any phase transition. This resembles the mass fraction trends 
within region 1 , only that the freeze-out within region 6 originates from equilibria states which are 
sensitive to the number and shape of clusters formed, and thus from the initial composition for the 
burning process. The physical border between regions 3 and 1 is the thin chasm line oriented ~ 70° 
with respect to the peak temperature axis. Such a distinctive border does not exist between regions 
1 and 6, due to the lack of a phase transition in both regions. However, an approximate border is 
the locus of points given by Tqse = 0.012 Tfreeze, where tqse is the timescale to reach QSE ( |Calder| 
et al.|2007 1 and Tfreeze is given by equation [3] This locus is shown by the thin cyan line in Figure |4j 



The relative differences for '*'*Ti yields starting from pure ^^Si or '^C are less than 0.1 to the right 
of this locus. 

The case of ^^Ni is simpler than '^'^Ti. The isotope ^^Ni tends to dominate the final composition 
for the majority of the peak temperatures and peak densities for 7^ = 0.5. The topology of its final 
mass fractions in Figure shows ^^Ni does not sustain any phase transitions like "^"^Ti because 



^^Ni remains in equilibrium with its local neighborhood (Woosley & Hoff'man 1992). While the 
macroscopic behavior of the large QSE cluster changes in different regions, there are almost no 
changes in ^^Ni. 

In region 1, a single QSE cluster which includes ^^Ni stays intact until freeze-out. In region 
2, the QSE cluster dissolves into two smaller ones, with the cluster localized around the Fe-group 
nuclei encompassing ^^Ni at all times. During the or-rich freeze-out of region 3 the QSE cluster 
shifts upwards in mass and shrinks ( Meyer et al.|1998 1, but remains centralized on Fe-group nuclei 



(including ^^Ni). Near the end of the evolution, the Fe-group nuclei are the most abundant in the 
network with reaction equilibria maintained among them. Figure [8] shows the mass fractions of 
^^Si, '^'^Ti and ^^Ni for a normal and an a-rich freeze-out, accompanied by the corresponding NSE 
values for each isotope were the NSE valid at all times. For the normal freeze-out, the network 
values are in good agreement with the corresponding NSE values until the NSE threshold of Tg ~5. 
For ^^Ni, the agreement between the network and NSE values persists until at least Tq ~2, at which 
point our NSE solver fails to converge. During an a-rich freeze-out the network values of ^^Si 
and "^"^Ti are quite different from their corresponding assumed NSE values, while the NSE mass 
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fraction of ^^Ni still agrees with reaction network values until ~2. Of course NSE at Tg ~2 
does not exist, but the trends in Figure [s] suggest that ^^Ni may be considered to be in a large scale 
equilibrium throughout the evolution for almost every region on the temperature-density plane. 
That is, global equilibrium estimates may interpret adequately the dominant trend of ^^Ni for an 
initially symmetric composition. 



5.1. Electron fraction sensitivity study 

The electron fraction, or the total proton to nucleon ratio. Ye = Y,i ^t^t = Z; f^^; » is equivalent 
to a weighted average of isotopic proton to nucleon ratios, where each has a probability equal to 
Xi. Since the distribution of isotopic ratios in a large network may be approximated as continuous, 
the most abundant isotopes at any given time in the thermodynamic evolution are generally the 
ones whose individual proton to nucleon ratio is within a small range from the current value of 
the electron fraction. A small spread usually exists due to nuclear structure effects for equilibrium 
states (expressed primarily with Q values), and reaction rate values for non-equilibrium states 
(Amett [T977| ). The largest (major) nuclear flows tend to be localized along the most abundant 



nuclei, since the flows depend on multiplications of abundances. During NSE or QSE the major 
flows result in the most robust reaction equilibria, while the same reactions typically become the 
most efficient carriers of nuclear flow as soon as they depart from equilibrium. Almost all our Ye 
sensitivity results may be explained by these guidelines for the major flows. An exception exists for 
cases with initial Yg > 0.5 during large scale equilibrium (NSE and QSE), where the equilibrium 
patterns are configured according to a different principle ( Seitenzahl et al.|2008[ ). Electron fraction 



variations alter the nuclear composition and affect the yields, the nucleosynthesis mechanisms for 
each region in the peak temperature-density plane, and change the regions topology. 



Figures M and 10 show the final yields of Ti and Ni, respectively, in the peak temperature 



density plane for 0.484 < Fg < 0.506 under the exponential freeze-out profile. Figure 1 1 shows 
the '^'^Ti production factor P44 for the same range of Ye- The production factor for a given species is 
defined as the final mass fraction of the species in question divided by the mass fraction to which 
it decays in the Sun. These production factors are then normalized to the production factor of ^^Fe 



( |Woosley & Hofi'manlligQlj [Hoffman et al.||2010| ). Within the electron fraction values 0.498 and 



0.5, yields for both isotopes are maximized, resulting to the minimization of the chasm's width 
and depth for '^'^Ti. For decreasing Ye values, both isotopes tend to be under-produced compared to 



the symmetric case (Woosley & Hofi'man||l992). For increasing Ye values, ^^Ni is still favored by 



equilibria schemes and is produced at an amount comparable to the symmetric case (Magkotsios 
et al. 2008). The temperature-density planes for ^^Ni have similar featureless structure to the 



corresponding plane for initially symmetric matter, implying that this isotope is produced only by 
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equilibria schemes without sustaining any phase transition. The location of the border between 
the regions of a-rich and orp-rich freeze-out for "^"^Ti depends on the initial electron fraction value. 
The lack of free protons for neutron-rich environments favors the a-rich freeze-out versus the 
ap-rich one, until the ap-rich freeze-out is not manifested at all for Ye « 0.46. The situation is 
gradually reverted for increasing 7^, until the ap-rich freeze-out dominates the a-rich freeze-out 
for Ye ~ 0.506. These trends are in accordance with the major flows guidelines discussed above, 
since both isotopes are symmetric with an individual proton to nucleon ratio equal to 0.5, and the 
amount of free protons increases significantly for Y^ > 0.5 ( Seitenzahl et al.|2008 ). Further changes 



to the topological structure compared to the symmetric case include the appearance of a depletion 
region for both "^"^Ti and ^^Ni for decreasing Ye, and the appearance of a physical border between 
regions 1 and 6 for increasing Ye. This physical border implies a new type of phase transition that 
"^■^Ti sustains. These trends are the same for the power-law freeze-out profile (not shown). 

We focus next on two relatively extreme Ye values, 0.48 for the neutron-rich case and 0.52 for 
the proton-rich one. Both values are adequately far from the standard value of symmetric matter, 
so that the differences in the trends for "^"^Ti and ^^Ni are emphasized and easily identified. The 
characteristic regions of the "^"^Ti temperature-density planes with initial electron fraction 7^=0.48 
and 7^=0.52 are labeled in Figure 12 for the exponential profile. Similarly to Figure |4} the only 



region sensitive to the initial composition is region 6, the incomplete Si-buming regime. 

For yg=0.48, NSE and QSE favor the formation of nuclei with a proton to nucleon ratio around 
0.48 during equilibrium, and the major flows are localized in the neighborhood of the same nuclei 
during the non-equilibrium parts of the evolution. Despite ^^Ni being marginally within the range 
of major flows, large scale equilibrium patterns gradually favor ^^Fe instead of ^^Ni for decreasing 
electron fraction (see Figure [T3|). Regions 1-6 each have the same type of physics compared to 
the corresponding ones for initially symmetric matter. The chasm widening is an outcome of the 
overall underproduction of '*"^Ti. The large scale equilibria patterns do not favor '^'^Ti production 



(see Figure 13), and the normal freeze-out region merges with the chasm. The a-rich freeze-out 
yields less "^"^Ti compared to the symmetric case, ceding additional area to the chasm region. The 
decreased efficiency of the a-rich freeze-out regarding '^'^Ti production is related to the manner that 
the electron fraction value affects the flow transfer by (a, y) reactions towards the N = 22 isotone 
chain, and the favor of the isotone chains towards neutron-rich isotopes rather than '^'^Ti. The size 
reduction of the ap-rich freeze-out region is due to the absence of free protons in neutron-rich 
equilibrium configurations. 

Region 7 represents the case of neutron-rich, a-rich freeze-out, which barely appears for 
Ye=0.5 (not labeled in Figure |4]). It is also known in the literature as the "a-process" ( |Woosley 



& Hoff'man 1992), but we term it henceforth as "an-rich freeze-out". This region combines the 



physics of regions 3 and 5. A photodisintegration regime is established early in the evolution and 
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during the equilibrium stages, but contrary to region 5, p(e", Vs)n and n(e^, y^p tend to balance each 
other after the phase transition imparted by the 3a forward rate dominance over its inverse. Thus, 
the electron fraction value is maintained close to its initial value, well below 0.5. Such values of the 
electron fraction are prerequisite for the production of elements beyond the Fe-group. Although 
the electron fraction has similar values within region 3, there are traces of Si-group and Fe-group 
nuclei during the equilibrium stages. The presence of ^^Ni blocks the flows towards heavier nuclei 
during the non-equilibrium stage and neutrons are consumed among the Si-group and Fe-group. 
However, these traces are absent for the equilibrium stages within region 7, and heavier elements 
are produced during the non-equilibrium stage. Table [2] lists the dominant yields from freeze-outs 
for Ye = 0.48, 0.50, and 0.52. The yields from the a-rich freeze-out are quite similar for Yg = 0.48 
and 0.50 ( Woosley & Hoffman|1992[ ). Since the region for the an-rich freeze-out increases in size 



for decreasing Y^, it is expected at some point to dominate all the area enclosed by the chasm. 
For 7^=0.52, NSE and large scale QSE are dominated by ^^Ni and free protons, with non 



negligible abundances for symmetric and proton-rich nuclei (also see Figure 13), in accordance 
with a minimum of the Helmholtz free energy. Despite the excess of free protons favored by 
large scale equilibria patterns, neutrons are captured more efficiently in a proton-rich environment. 
Thus, n(e"^,y^)p reaction dominates p(e~,Ve)n, resulting in a slight increase to the value of the 
electron fraction early in the evolution for the cp-rich freeze-out region. Since these reaction 
channels retain Y^ > 0.5 during large scale equilibrium, major flows favor proton-rich nuclei as 
soon as the QSE cluster dissolves. Thus, for regions 1-5 where NSE and large scale QSE are 
established, the (p, y) reactions are all directed towards transferring the flow from symmetric to 
proton-rich isotopes, which is equivalent to a phase transition that all isotopes in the network 
sustain. Consequently, there cannot be a normal freeze-out in region 1 (Figure [12]), because the 
freeze-out does not take place from NSE (or large scale QSE). During the non-equilibrium part 
of the freeze-out evolution, the weak interactions decrease Y^ by transferring flow towards more 
stable isotopes, resulting in its non-monotonic evolution and the reassemble of symmetric isotopes 
hke 44Ti. 



Within the ap-rich freeze-out region (region 4 in Figure 12) the '^'^Ti mass fraction pattern 
resembles the corresponding one with initial Y^, = 0.5, with two arcs and an ascending track at the 
end. A timescale dependent third arc is identifiable for the power-law profile only. Its appearance 
relies on the equilibrium state of the remaining '^^Ti-'^^V cluster and the net flow towards this 
cluster by the interplay between neighboring (p, y) and weak reactions. V(p, y)'^^Cr is the primary 
reaction to control the flow leakage off this cluster. 

Within region 6, the initially formed small scale QSE clusters fail to merge to a large-scale 
cluster. This results in randomly directed flow supply among the small scale clusters and the 
absence of a phase transition accompanied by complete consumption of fuel nuclei. The physical 
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border between regions 1 and 6 is an outcome of the existence of a transition within region 1 . For 
regions 1-5 the final composition is always proton-rich. 



6. Reaction rate sensitivities 

The topology of the "^"^Ti and ^^Ni contour plots is affected by certain key reactions, in com- 
bination with the timescale of the expansion. We follow a three-stage method to uncover the role 



and impact of these reactions. Figure 14 exemplifies the three stages of this method. During the 
first stage, specific reaction channels are either altered or removed from the network calculations 
for all isotopes (e.g., all (a, y) reactions) to assess the most significant channels for every region. 
In addition, the 3a, p(e", Ve)n and n(e^, v^)p reactions have their own brevet, since they may affect 



the reaction flows globally (first row in Figure 14). Thus, the term "weak reactions" will imply all 
such reactions henceforth, excluding p(e , Ve)n and n(e^, v^)p. We tabulate weak reactions by their 
dominant decay mode, although all decay modes are considered in our calculations. The second 
stage performs a sensitivity analysis on all groups of reactions by increasing and decreasing reac- 



tion rates excessively one at a time (third panel within first row in Figure [14]). Similarly to |The 



et al. ( 1998), reaction rates are either multiplied or divided by a factor of 100. The exception are 
the weak reactions, where the factor is 1000. The third stage conducts detailed nuclear flows and 
mass fraction profile analysis to illustrate the impact of the final crucial reactions that affect "^"^Ti 



(second row in Figure 14). Note the second and third stages are applied independently for every 



distinctive thermodynamic region in the temperature-density plane. The rates used for our calcula- 



tions are from the ,Rauscher & Thielemann] ( jlOOOj ) compilation, updated with some experimentally 



measured rates. Table[3]lists the most important reactions which our sensitivity study has revealed 
to impact '^'^Ti. We rank reactions as "primary" or "secondary", depending on the differences be- 
tween the "^"^Ti mass fraction curves for nominal and modified rates. A reaction which involves 



differences at any point of the evolution by a factor of 10 or larger is ranked primary (Figure 14). 
Reactions resulting in changes smaller than a factor of 10 are ranked secondary. Reactions of 
minimal impact are not tabulated. 

It is important to clarify the advantages and disadvantages of our methodology for the sensi- 
tivity studies. Within the regime of medium mass nuclei where "^"^Ti and ^^Ni belong, the nuclear 
level densities are large enough, so that uncertainties to reaction rates are expected to be constrained 
within a small range from their nominal values ( Iliadis|2007 ). However, such small changes to the 



rates may not fully demonstrate the impact of individual reactions to the burning process. Our 
goal is to understand the microscopic mechanisms of explosive nucleosynthesis, which are driven 
by the effect of individual reactions in combination with localized equilibria patterns. Unrealistic 
changes to reaction rates either by excessive factors or by removal from the network are required 
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to result in distinguishable changes to the dynamics of the burning process. The changes to the 
burning process are related to isolated microscopic components to the operation of explosive nu- 
cleosynthesis. Our sensitivity study aims to identifying as many as possible of the components 
related to "^"^Ti and ^^Ni synthesis. Thus, we add detail to our understanding of the process for nom- 
inal values of the reaction rates. On the other hand, this type of sensitivity study may not provide 
a numeric measure of the importance of identified reactions. For this purpose, sensitivity studies 
should be constrained within acceptable uncertainty limits for the reaction rates. Such sensitivity 
studies have been performed by [Hoffman et al.|(|2010[) and|Tur et al.|(|2010l). 



Overall, the (p,n), (a,n) and (n,7) reactions have either secondary or minimal impact to the 
synthesis of '*'*Ti. The reason is that '^'^Ti is produced mostly for Yg > 0.5, where neutrons tend to 
be depleted quite fast. Moreover, reactions that emit neutrons usually have higher thresholds than 
proton emitting ones, since neutron separation energies are larger than proton separations energies 
for proton-rich nuclei. 



6.1. The 3 a reaction 

One difference between the normal and a-rich freeze-outs is the behavior of the large scale 
QSE cluster. For decreasing temperatures and densities local equilibria successively break, gradu- 
ally dissolving the cluster. The most sensitive equilibria are related to reactions with large Q values. 
When no external flows are applied to the QSE cluster, the first equilibria to break are among iso- 
topes with 42 < A < 53, where the largest reaction Q values for alpha particle captures in the 
network are localized (Figure [T]). This is the case for the chasm, region 2 in Figure |4j When these 
local equilibria break, the QSE cluster dissolves into two smaller QSE clusters; the first localized 
within the Si-group nuclei, and the second within the Fe-group nuclei. During this process, the 
thermodynamic conditions dictate the net 3a rate is always dominated by its photodisintegration 
reverse rate. 

In contrast, the forward flow dominates the net 3a rate for thermodynamic conditions con- 
ducive to an ff-rich or ap-rich freeze-out. Here the 3a reaction supplies the external flow to the 
large equilibrium cluster from the region of light nuclei. Specifically, reactions in the neighbor- 
hood of ^"^Mg have relatively large Q values (Figure [T]), although slightly lower compared to the 
ones in the region 42 < A < 53, and are the first equilibria to break under contributions from the 
3a reaction. The external flow supply results in a phase transition, leading to the a-rich freeze-out. 
Omission of the 3a reaction from the network calculations results in the severe underproduction 



of the Si- group elements as shown in Figure 14 for Ti. This happens because the phase transi 



tion is prohibited from taking place, and freeze-out from QSE at these conditions favors only the 
Fe-group nuclei. Omission of the 3a reaction within the normal freeze-out regime has little effect 
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as the forward rate has no impact in this regime. 



6.2. The {a,y) reactions 

Following the phase transition during an or-rich freeze-out, alpha captures break equilibrium 
and transfer nuclear flow between (p, 7) equilibria chains along isotone lines. Depending primarily 
on (i) mass differences between reactants and products of a reaction and (i) the electron fraction 
value, (or, 7) and (a, p) channels compete for the dominance in flow transfer. The Q value for 
'^*'Ca(a,p)'^^Sc (Q ~ -3.522 MeV) allows the gradual dominance of "^^CaCo-, 7)'^'*Ti for decreasing 
temperature in both the a-rich and orp-rich freeze-out regions. The impact of this reaction appears 
as soon as "^"^Ti moves off the large scale QSE cluster, after the first dip in the mass fraction curves 
caused by (a, p) and (a, 7) equilibrium breakages (Figure [T?]). 



In accordance with the major flow guidelines, '^*'Ca(Q', 7)^'*Ti is the primary reaction to supply 
flow along the N=22 isotone for symmetric matter. This supply is responsible for maintaining the 
pattern of (p, 7) equilibria along that chain for decreasing conditions. Were this flow supply absent, 
the (p, 7) equilibria chain would break and the ascending track of "^"^Ti mass fraction would cease. 
Thus, '^^Ca(a, 7)'*'^Ti regulates the amplitude of the subsequent rise past the first dip. Breakage 
of various (p, 7) equilibria determines the formation of additional such dips, and '^"Ca(a, 7)'^'*Ti 
regulates the amplitude of the formed arc in the mass fraction curve. For the a-rich freeze-out of 



Figure 14, "^^Ca^a, 7)'*'*Ti controls how high the mass fraction curve rises once past the dip. Larger 
rates enhance the flow into the N=22 isotone chain, resulting in an increase in the '^'^Ti yield. 



However, a larger rate has the opposite effect in the ap-hch freeze-out of Figure 14 The second 
phase transition is caused by "^^VCp, 7)'^^Cr breaking from equilibrium (see § 6.4 1. The depth of the 



second dip and the magnitude of the subsequent ascent is controlled by '^^Ca(a, 7)'^"^Ti. A larger 
'**'Ca(a, 7)'*'*Ti rate enhances the depth of second minimum, resulting in a smaller overall "^"^Ti yield. 



Figure 14 shows that '^°Ca(Qr, 7)'*'*Ti affects the amplitude of the second arc for the '^Ti mass 
fraction, but the slopes of the ascending and descending tracks are relatively robust. These slopes 
are determined by the (p, 7) channels that participate in the equilibrium chain. Their robustness for 
symmetric matter is a direct consequence of the major flows guidelines. However, proton captures 
are less efficient within a neutron-rich environment, and the slopes are affected by the net flow 
transferred to the mildly connected equilibrium chain. The net transfer is determined primarily 
by the flow supply from '**'Ca(a, 7)'^'^Ti and '*°Ca(ar, p)'^^Sc and the flow leakage from '^'^Ti(a, p)'^^V 



(see also ^6.3) 



For Ye < 0.5 the impact of '*'^Ti(a, p)'*^V on the chasm (see §6.3) is influenced by the sec- 
ondary '^^Ca(a, 7)'*''Ti. Further minimal contributions from other {a,y) reactions are related to the 
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distribution of nuclear flow among the remaining equilibria chains along various isotone lines past 
the QSE cluster dissolution. In addition, '^C(a, 7)^^0 is a secondary reaction to affect the flow 
supply to the QSE cluster after the 3a rate has dominated its inverse, with minimal contributions 
from 20Ne(a, r)^^Mg. 

For Ye > 0.5, the amplitude regulation of the second arc by ^°Ca(Qr, 7)'^'^Ti has an impact to 
the "^^Ti yield only for the exponential profile, due to the absence of a third arc in the "^"^Ti mass 
fraction for this case. The rest of the (a, y) reactions are secondary to "^"^Ti synthesis for initially 
proton-rich composition. 



6.3. The (Qr,p) reactions 

Of vital importance to "^^Ti synthesis from this channel group is '*'^Ti(a', p)^' V. This reaction 
is related directly to the formation of the chasm, which is the border region between the normal 
and the a-rich freeze-outs (Figure |4]), the depth of the chasm, and the location of the chasm in 
the peak temperature-density plane for different expansion profiles. However, this reaction is not 



responsible for the widening of the chasm; weak reactions discussed in §6.5 largely control the 
chasm width. The key feature of '^'^Ti(a,p)'^'^Y is its small negative Q value (Q ^ -410 keV). 
This feature allows '*'*Ti(Q', p)'^'V to dominate '^'^Ti(Q', 7)'^^Cr even for low temperatures. Based on 
the major flow guidelines, '^'^Ti(a, p)'*^V is the primary flow supplier from the N=22 to the N=24 
isotone for initially symmetric matter. When this reaction is in equilibrium, '^'^Ti is considered to 
belong in the large scale QSE cluster, a fact verified by QSE calculations (Figures [6] and |7]). Its 
equilibrium breakage signals the phase transition for '^'^Ti, leaving the isotope outside the QSE 
cluster. 

The '^'^Ti chasm is formed by dissolution of the large QSE cluster into two smaller clusters, 
and the subsequent flow leakage from one cluster to another. For peak temperatures and densities 
corresponding to the chasm region, '*'^Ti(Q', p)'^'V is always in equilibrium until the very end of 
freeze-out. The resulting mass fraction for "^^^Ti ends up with a yield 2-3 orders of magnitude less 
than its typical value in regions outside the chasm region, as shown in Figure |5j For the or-rich 
and ap-rich freeze-out regions, '*'^Ti(Q', p)^^V breaks equilibrium before the end of the freeze-out. 
Equilibria transitions for '^'^Ti are depicted in Figure |6] There is a robust equilibrium between the 
two isotone chains initially, but eventually '*'*Ti(a, p)"^^ V departs from equilibrium. Although this is 
an endothermic reaction, the a capture dominates its inverse because free alpha particles are more 
abundant than free protons {X(a) » X(p)). The rest of the equilibria connecting the N=22 and 
N=24 isotone chains break sequentially. When no equilibria links connect the isotone chains, the 
abundances of all related elements begin to increase. From this perspective, the '^'^Ti(a, p)'^^V reac- 
tion's persistence until the end of freeze-out is important for all isotopes along the N=22 isotone 
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chain, not just "^"^Ti. 

The chasm's depth is directly related to the minimum value of the mass fraction curve for "^"^Ti 
prior to the equilibrium breakage of ^^Ti{a, p)'^' V. Sensitivity studies for this reaction reveal that 



the minimum value is determined by the rate's strength. In Figure 14 the "^^Ti mass fraction during 
an a-rich freeze-out is shown as a function of temperature for various multiplicative factors to the 
'*'^Ti(a;, p)'^^V rate. The minimum value is smaller for larger reaction rates. Note that the slope of 
the mass fraction curve after the minimum value is independent of the rate's strength, showing this 
reaction has no impact on '^'^Ti synthesis from the moment this reaction goes off equilibrium. 

One of the major differences between the exponential and power-law profiles is the location 
of the chasm in the peak temperature-density plane. Figure [3] shows the chasm occurs at smaller 
densities for the power-law profile. The reactions that change the yield of '^'^Ti between these two 
profiles are approximately the same. This excludes reactions alone as a reason for the location of 
the chasm, implying timescale effects play a key role. The power-law expansion always evolves 
slower than the exponential one for the same initial peak temperature and peak density. In an 
environment where nucleosynthesis is driven by entropy changes, temperature sets to first order 
the threshold for a particular phase transition to appear, but the density value at the threshold 
temperature determines whether the transition takes place or not. The contribution of timescale 
effects to the chasm shift is related to the time spent by the plasma in between phase transitions 
(Figure [2]), resulting in different density values at threshold temperatures. Thus, the chasm shift is 
a density driven phenomenon. Specifically, '^'*Ti(Q', p)'^'V departs from equilibrium approximately 
at the same temperature Tthr ~ 4.3 GK for both expansion profiles. When both the exponential 
and power-law profiles reach Tthr, the density associated with the exponential profile is larger (and 
earlier in time) than the density of the corresponding power-law profile (later in time). Since both 
profiles assume a constant radiation entropy, p « T^, throughout the evolution, a larger (or smaller) 
density at Tthr translates directly into a larger (or smaller) initial peak density. This causes the shift 
in the location of the chasm in the peak temperature-density plane. 

For neutron-rich environments, '^^Ca(a, p)'^^Sc pipes flow from the major flows among neutron- 
rich isotopes to "^"^Ti, increasing thus its mass fraction. The reaction's main feature is a flow di- 
rection switch for conditions past the phase transition. While the "^"^Ti mass fraction begins its 
ascending track, the forward flow of '^*'Ca(a, p)'^^Sc dominates its inverse. When the flow direction 
for '*°Ca(ar, p)'^^Sc switches and the proton capture dominates the alpha capture, part of the ma- 
jor flows is supplied to "^"Ca, and subsequently to the N = 22 isotone equilibrium chain through 
'**'Ca(cir, y)'^'^Ti. This pattern of escalating flow exchange between the N = 20 and N = 22 isotone 
chains has also an impact within a proton-rich environment. During the formation of the second 
arc for the '^'^Ti mass fraction, the dominant proton capture in '^*'Ca(a, p)'^^Sc results in an enhanced 
flow supply to '^'^Ti, which is lost during the leakage to proton-rich nuclei by '^^ V(p, 7)'^^Cr. Hence, 
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the Ti mass fraction decreases in this case. 



6.4. The (p,7) reactions 

This group of channels is characteristic for the collective contribution of reactions in the form 
of equilibria chains. The most important proton captures are localized among symmetric and 
proton-rich isotopes, because their inherently enhanced efficiency may alter equilibria patterns 
and result in phase transitions. Their effectiveness is enhanced significantly in proton-rich envi- 
ronments, where they are favored by major flows and there is a large availability of free protons. In 
practice, the weak reactions set a proton-rich environment primarily with p(e", Ve)n and n(e^, v^)p. 
Without this elegant combination of weak interactions and proton captures, the ap-rich freeze-out 
region in the contour plots merges smoothly with the a-rich freeze-out one. Proton captures al- 
ter the local equilibria patterns, resulting in small scale phase transitions. For initially symmetric 
matter, they sculpt the ap-rich freeze-out topology in the contour plots. 

The (p, y) channels most relevant to '^'^Ti nucleosynthesis operate along the A'^ = 20, N = 22 
and N = 24 isotone chains. The important isotone chain is the N = 22 one, where '^'^Ti resides. 
The specific chain is composed by '^'^Ti, "^^V, '^^Cr, "^^Mn and "^^^Fe. The chain terminates to "^"^Ti 
due to the early equilibrium break of "^^ScCp, 'y)'*'^Ti, while the upper limit of '^^Fe appears due to 
the large negative Q value of '^^Fe{p, yy^Co close to the proton dripline. For regions 3 and 4 in the 
temperature-density planes, the ascending part of the second arc in the mass fraction profile for '^'^Ti 
is formed when these isotopes are all in mutual equilibrium. For > 0.5, the major flows attribute 
a relative robustness to the slope of the ascending track from single reaction sensitivities. This 
robustness is gradually fading as material becomes neutron-rich, and the equilibrium maintenance 
along the chain depends on the net flow supply, which is configured primarily by '^"Ca(Q', y)'^'*Ti, 
^'^Csiia, p)43Sc and ^'^Ti{a, p)^'V. 

Among the reactions connecting the isotopes within the N = 22 isotone equilibrium chain, 
"^^VCp, y^^Cr has the largest Q value, rendering it the most sensitive equilibrium link. Within the 
ap-rich freeze-out region, it is the first one to break, leaving '^'^Ti in equilibrium only with ''^^V and 
prognosticating a phase transition where the mass fraction of "^"^Ti decreases ([The et al.|l998 1, due 



to the flow transfer from the '^'^Ti-'^^V cluster to the '^^Cr-'^^Mn-'^^Fe cluster. At the same time, the 
equilibria patterns along the rest of the related isotone chains change, contributing all together to 



the phase transition for '*'*Ti. Furthermore, '^^V(p,y)'^^Cr is secondary (see §6.5 below) to defining 
the physical border between the regions of a-rich and ap-rich freeze-outs. A stronger rate expands 
the ap-rich freeze-out region at the loss of the or-rich freeze-out region. The ^^Ni(p, 7)^^Cu is 
another secondary reaction which contributes to the localization of the physical border between 
regions 3 and 4. When the large scale QSE cluster begins to dissolve, it is one of the primary 
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reactions to control the flow transfer within the remnant QSE cluster. 

A couple of reactions with a sensible impact to the "^"^Ti yield are "^'ScCp, y)^^Ti and '^^ScCp, 'y)'^^Ti. 
They are the primary reactions to regulate the depth of the second dip in the "^"^Ti mass fraction, 
by affecting the flow supply to the equilibrium chain by ^Ca{a, y)'^'^Ti. The '^"^Ti(p, y)'^^V reaction 
is the immediate link of "^"^Ti to the specific equilibrium chain. A stronger rate maintains the ex- 
istence of the '^'^Ti-'^^V cluster, resulting in further loss of flow via '^^V(p, yy^Cr. Thus, the "^"^Ti 
yield is decreased. A secondary reaction to affect the ascending track beyond the second arc is 
^''Ca(p,r)4iSc. 

For proton-rich environments, the (p, y) channels are primary to the formation of the (p, y)- 
leakage region (Figure [12]). Large scale equilibria patterns favor both symmetric and proton-rich 
nuclei ( Seitenzahl et al.|[2008| ). As soon as the large scale QSE cluster begins to dissolve, (p,7) 



reactions transfer the flow from symmetric nuclei to proton-rich ones, so that the major flows are 
localized in the neighborhood of the latter, in accordance with the major flows guidelines. Without 
the contribution of the (p, y) reactions, the (p, y)-leakage region would be equivalent to a normal 
freeze-out regime and would merge smoothly with region 6, such as the = 0.5 and Yg < 0.5 
cases. The flow transfer by (p,y) reactions is massive, where almost all of them in the network 
participate. Thus, the initially descending track of the '^'^Ti mass fraction due to this flow transfer 
is relatively robust to single rate sensitivities. The '*^V(p,y)'*^Cr reaction is the only one to affect 
the depth of the descending track, especially for the power-law expansion profile. 

In addition, "^^V(p, y)'*^Cr controls the flow leakage off the remaining ^^Ti-'^^Y cluster during 
the ap-rich freeze-out, once this reaction breaks equilibrium. In combination with the timescale 
of certain weak reactions in the locality of "^^^Ti this results in the formation of the third arc for the 
"^"^Ti mass fraction for the power-law profile. Secondary reactions within this group (along with the 
weak reactions) which regulate the amplitude of the third arc are listed in Table |3j 



6.5. Thep(e , Ve)n and n(e^, Ve)p and weak interactions 



The electron fraction Ye expresses the proton to baryon ratio in the plasma. Assuming charge 
neutrality, the electron fraction is also the electron per baryon ratio. Weak interactions are the only 
group of channels to violate the lepton number conservation, while preserving the baryon number. 
Thus, they are the only ones to change Y^. during the evolution, with p(e~, Ve)n and n(e"^, v^)p having 



a special contribution to this configuration (McLaughlin & Fuller 1995; 


Fuller & Meyerl 1995; 


McLaughlin et al.|1996 


Surman & McLaughlin|2005 , Liebendorfer et al.|2008[ Aprahamian et al. 



2005]). Depending on the competition between p(e ,Ve)n and n(e"^,ye)p, the electron fraction Yg 
may increase or decrease. For these two reactions to be effective, relatively large temperature and 
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density values are needed. Thus, their impact is usually constrained during the first stages of the 
evolution. On the contrary, the lifetimes of the remaining weak interactions ensure their impact 
appears during the last stages of the evolution. These reactions tend to transfer material towards 
the valley of stability. Our calculations use the FFN rates for the p(e~,Ve)n, n(e"^,y^)p and other 
weak reactions ( jFuUer et al.|1980[ |1982b|a| |Oda et al.|[T994] [Langanke & Martmez-Pinedo|2001| ). 



Using the FFN weak rates for the other reactions has little effect on the synthesis and yields of ^"^Ti 
and ^^Ni. We thus use temperature and density independent jS -decay and ;S"^-decay rates, where 
the parent nucleus is assumed to be on its ground state. 



In combination with the nucleosynthesis trends for a varying electron fraction (see § |5.1| 
above) and timescale effects, the p(e",Ve)n and n(e^,v^)p reactions are the key to explaining 
the chasm's widening between the exponential and power-law expansions. Depending on the ex- 
pansion timescale, p(e", Ve)n and n(e^, v^)p reactions alter the electron fraction only for a limited 
amount of time early in the evolution. The changes to for various peak temperatures and densi- 



ties depend on the rate strengths of these reactions and affect the yields directly. Figure 15 shows 
the evolution of and a few isotopes related to "^^Ti nucleosynthesis for a case of a normal freeze- 
out from initially symmetric matter (Y^ = 0.5) using nominal rates for both expansion profiles. In 
this regime, p(e^, Ve)n always dominates n(e^, v^)p and Ye decreases, while temperature and density 
still have large values. Timescale effects are evident when using nominal rate values, where the 
time spent in a high entropy environment is larger for the power-law case and Y^ decreases much 
more compared to the exponential profile. The plasma adjusts to the Y^ 4^ 0.5 conditions while it 



is still in NSE and QSE subsequently. Figure 13 shows NSE mass fractions as a function of the 



electron fraction where the production of ^^Fe is favored, while "^^Ti and ^^Ni are under-produced 
Partmann et al.||1985||Woosley & Hoffman|1992t|Seitenzahl et al.||2008D . On the contrary, expo- 



nential expansion does not allow Y^ to decrease significantly, resulting in a final composition with 
significant yields for "^"^Ti and ^^Ni. Therefore, the chasm expands only for the power-law profile. 

The chasm width is regulated primarily by the strength of p(e^, Ve)n and n(e^, y^)p and secon- 
darily by timescale effects. Both '*'*Ti and ^^Ni show a large chasm expansion for both thermody- 
namic profiles when these two reactions are enhanced by a factor 10^. The new chasm widths for 
the two profiles are similar, because the reaction rate enhancement results in the same decrement 



to Ye for both thermodynamic profiles and diminishes the impact of timescale effects. In Figure 14 
the normal freeze-out regime (region 1) for "^"^Ti merges with the chasm (region 2) and the chasm 
expands into the area that belonged to the or-rich freeze-out regime (region 3). 

The remaining weak interactions also assist in the decrement of the electron fraction, and thus 
to the chasm expansion, but their contributions are smaller than p(e , Ve)n and n(e^,y^)p due to 
their lifetimes. The lifetime of any weak interactions that are primarily responsible for the changes 
to the electron fraction must be smaller than the expansion timescale. In Figure [15] the changes to 
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Ye take place between 10"^ < t < 10"^ sec. The exponential profile has a timescale of the order 
of 1 sec, and the changes to Ye are moderate. On the contrary, the timescale for the power-law 
is larger by two orders of magnitude, resulting in dramatic changes to Ye due to the impact of 
p(e~, Ve)n, n(e"^, v^)p, and weak interactions. Despite the initial identical configuration of the two 
expansions, the equilibrium state that the normal freeze-out begins is very different for the two 
expansion profiles. For the exponential trajectory, the weak interactions do not have the time to 
change the equilibrium state adequately, and the final yields have a significant amount of '^Ti, with 
^^Ni dominating the final composition. For the power-law profile ^^Fe is the dominant element and 
"^"^Ti is under-produced, expanding the chasm region into the normal freeze-out regime (see Figure 



15[ ). In addition to p(e ,Ve)n and n(e^,ye)p reactions. Table |3J lists the primary weak interactions 
related to the chasm widening, all with a half life of the order of 10"' sec. 

Weak interactions assist in defining the '^''Ti topology for the ap-rich freeze-out regime (re- 
gion 4). In this region, n(e'^,v^)p dominates over p(e",Ve)n and Ye rises above 0.5, driving the 
material proton-rich. The relative strength of the p(e , Ve)n and n(e'^,v^)p rates determines the 
area of the peak temperature-density plane occupied by the ap-rich freeze-out regime as shown 



by Figure 14 Both '^'^Ti and ^^Ni have significant mass fraction values during equilibrium states 
in such environments, despite the relatively large mass fractions of proton-rich isotopes. These 
isotopes decay within the expansion timescale and transfer additional nuclear flow to the symmet- 
ric isotopes. Weak reactions partially regulate the second arc for the '*'*Ti mass fraction, and the 
formation of the ascending track at the end of the freeze-out process. 

For Ye = 0.48 the weak interactions barely have an impact on the "^''Ti yield. Within the 
ap-rich freeze-out regime their action is similar to the symmetric case, but the area that region 
4 occupies on the temperature-density plane for neutron-rich matter is limited. Within the or-rich 
freeze-out regime the major flows are localized mostly among stable nuclei, or nuclei with decay 
timescales much longer than the expansion timescale. 

For Ye = 0.52 the action of the weak interactions has been outlined before. They trans- 
fer nuclear flow from proton-rich to symmetric nuclei for most of the peak conditions within the 
temperature-density plane. For the ap-rich freeze-out region with the exponential profile, the av- 
erage half-life range for the primary flow carriers is 90 < ?i/2 < 500 ms. For the power-law 
expansion, the corresponding range is 200 < ?i/2 < 900 ms. Weak reactions for the (p, y) leakage 
regime (region 1 for Ye = 0.52) are classified according to the way they impact. The first group 
includes reactions which hinder the flow transfer by (p, y) reactions when their rates are enhanced. 
The second group includes reactions which boost the flow transfer by (p, y) reactions when their 
rates are diminished. The third group includes reactions which combine the action from the two 
previous groups. Reactions within the third group make an impact only for the power-law expan- 
sion profile. In addition, weak reactions with relatively long half-hves contribute to the amplitude 
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regulation of the third arc in the '^'^Ti mass fraction for the ap-rich freeze-out regime (see Table [b]). 



7. Network size effects 

Trends in the "^"^Ti yields are controlled by a limited number of reactions. This raises a query 
about the minimum number of nuclei that are necessary to include in a network calculation such 
that all relevant physical phenomena are captured. Our reference reaction network contains 489 
isotopes, spanning the light nuclei, silicon group, and iron group. To assess network size ef- 
fects, we compared the '^'^Ti yields in the peak temperature-density plane from the 489 isotope 
network with the final yields generated by reaction networks with 204, 1341, and 3304 isotopes. 
Table [T] lists the isotopes used in each network. In addition to the new elements introduced, the 
larger networks expand into larger both neutron-rich and proton-rich regimes. The addition of 
new elements beyond the Fe-group has a minimal effect on equilibrium clusters, since for the 
thermodynamic trajectories of interest the largest partial flows are localized around the Si-group 
and Fe-group. 

For Ye = 0.48, the final mass fractions are essentially independent of the network size since the 
major flows occur along the valley of stability, which is modeled adequately by all networks. For 
Ye > 0.5, the final mass fractions depend on the network size as weak interactions have a larger role. 
In particular, the 204 isotope network does not include most of the required proton-rich isotopes 
to accurately describe the nucleosynthesis. The differences compared to our reference network are 
localized to the ap-rich freeze-out region for Yg = 0.5, but they span all the parameter space for 
Yg = 0.52. Consequently, this 204 isotope network is inadequate to describe the nucleosynthesis 
of proton-rich material. 

There are mass fraction differences in region 1 for Y^ > 0.5 between our reference network 
and the larger networks, which is related to the differences in the equilibrium state configurations 
by the changes in temperature, density and Ye during the freeze-out evolution. Specifically, the dif- 
ferences are due to the relationship between the expansion and weak interaction timescales. Figure 



16 shows the temperature dependence of the electron fraction during a freeze-out for the 489, 
1341, and 3304 isotope networks for the exponential and power-law trajectories. The 489 and 1341 
isotope networks have relatively similar numbers of isotopes per element, for those elements that 
are common to both networks. Consequently, the evolution of Ye is quite similar for both profiles. 
The 3304 isotope network, however, has a larger number of isotopes per element for elements that 
are common among the three networks (see Table [T]). The presence of more proton-rich nuclei in 
the 3304 isotope network slows the electron fraction decrease driven by p(e", Ve)n and n(e^,y^)p 
compared to our reference network. This results in slightly different large scale equilibrium states. 
For a short expansion timescale, such as the exponential profile. Ye values remain above 0.5 for all 
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networks, resulting in the differences in region 1. The long timescale of the power-law expansion 
decreases the electron fraction value below 0.5 quite early in the evolution. This results in all three 
networks converging to the same Yg values since all three networks include the necessary isotopes 
related to production of '^"^Ti and ^^Ni. Overall, our reference 489 isotope network is adequate 
for describing the trends in the "^"^Ti and ^^Ni yield trends. This is relevant for efficient use of 
computational resources. 



8. Post-Process Yields from Collapse Simulations 

In this section we compare the '*^Ti and ^^Ni yields from post-processing core-collapse su- 
pemovae models with the exponential and power-law trajectories. We use the same reference 489 
isotope reaction network for the post-processing and parameterized trajectories. Our aim is to offer 
a calibration of where parameterized trajectories provide a reasonable approximation to the final 
yields. Our analysis in the preceding sections allows an explanation for the behavior of the "^"^Ti 
and ^^Ni profiles and any differences between the post-processed and parameterized yields. For 
this assessment we consider 3 of the supernova explosion calculations whose tracks in the peak 
temperature-density plane are shown in Figure |3j In all three of these models the initial profile 
as a function of interior mass is very close to Ye = 0.5. 



8.1. A Cassioppeia A model 

Our first supernova model uses a progenitor designed to match the supernova remnant Cas 
sioppeia A ( Young et al.|2006 , Eriksen et al.|2009] ), specifically model M16El.lBinA from Young 



et al.| ( |2008| ). The hydrogen envelope of this spherically symmetric 16 M© progenitor was removed 
by an assumed binary mass transfer event as the progenitor evolved into a giant. We use a multi- 
step collapse and explosion process to model the explosion (for example, see Young & Fryer|2007 1. 



We model the entire star from collapse through the formation and stall of the bounce shock. At this 
point, the proto-neutron star is removed from the simulation and energy is injected just above the 
proto-neutron star to drive an explosion. The explosion is followed as the shock moves out of the 
star and most of the fallback has accreted onto the newly-formed neutron star. We only calculate 
the yields of material that is ejected after fallback. 

Mass fraction profiles from post-processing the Lagrangian thermodynamic trajectories with 
our reference 489 isotope network are shown in the top panel of Figure[l7] The iron-group, silicon- 
group, oxygen-rich shells are visible within the innermost 1.5 Mq. The bottom panel shows the 
mass fraction profiles from the post-process, exponential, and power-law trajectories. Mass frac- 
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tion profiles correspond to the left-hand y-axis, and the peak temperature and peak density curves 
correspond to the right-hand y-axis. Figure[3]also shows the peak temperatures and densities of this 
explosion model, and explains the general trends in the '^'^Ti, ^^Ni and "^He mass fraction profiles of 



Figure 17 



For interior masses less than « 0.2 Mq, the synthesis of '^'^Ti is due to the or-rich freeze-out, 
region 3 in Figure |4] In this mass range, the "^"^Ti mass fractions from the power-law profile are 
closer to the post-process values than the exponential profile. The ^^Ni mass fractions given by 
the power-law and exponential profiles generally agree with the post-process mass fractions. Near 
0.2 Mq, the thin chasm region separating region 3 (a-rich freeze-out) and region 6 (silicon-rich) 
is traversed, which causes the downward spike in the '^'^Ti profile (see Figures [3] and |4]). Precisely 
where the thin chasm is traversed depends sensitively on the exact values of peak temperatures 
and densities, and the location of the thin chasm line on the temperature-density plane due to the 
timescale of the expansion profile. This explains why the downward spike occurs at slightly dif- 
ferent mass locations for the post-process, exponential and power-law profiles. Between « 0.2 M© 
and X 0.5 Mg the explosion is operating in the silicon-rich, region 6 of Figure |3| In this region 
the final yields of both '^'^Ti and ^^Ni are sensitive to the initial composition, which in this case is 
given by the model at the time when energy is injected just above the proto-neutron star to drive an 
explosion. In this mass range, the '^'^Ti mass fractions post-processing, exponential, and power-law 
profiles generally agree. For ^^Ni, the power-law profile generally agrees with the post-process 
mass fractions better than the exponential profile. Both "^"^Ti and ^^Ni abruptly decline at « 0.5 Mq 
as the peak temperature drops below 4x10^ K. 

Yields from the exponential and power-law trajectories are generally within a factor ~ 2 of the 
post-process yields, except in region where the thin chasm is being crossed or the temperature falls 



below 4x10^ K. Integrating the "^Ti mass fraction profiles in Figure 17 over the interior mass gives 
the total mass of '^'^Ti ejected by this model. We find 1.04x10"^ Mq for post-processing, 5.62x10"^ 
Mq for the exponential profile, and 9.30x10"^ Mq for the power-law profile. Similarly for ^^Ni, we 
find 2.46x10'' Mq for post-processing, 3.16x10"' Mq for the exponential profile, and 2.78x10"' 
Mq for the power-law profile. Overall, yields of ''''Ti and ^^Ni from the power-law profile mimic 
the post-process values better than the exponential profile for this Cas A model. 



8.2. A Weak-Strong Hypernova Model 

Our second massive star explosion model uses a similar multi-step process. In this case how- 
ever, first a weak explosion is launched and followed 1 s later by a strong 1.6 x lO^^erg hypernova 
explosion. This model has a very different thermodynamic evolution as the material is hit by two 
shocks. The first, weaker shock ignites a substantial amount of burning within the inner 0.6 Mq. 
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Most of the peak temperatures and peak densities inside 0.6 Mg are due to the weaker shock. Peak 
conditions at larger masses are due to the second, stronger shock. For additional details on this 
weak-strong double shock model, see 40WS1.0 from Fryer et al. ( |2006| ). 



Mass fraction profiles from post-processing the Lagrangian thermodynamic trajectories with 



our reference 489 isotope network are shown in the top panel of Figure 18 The iron-group and 
silicon-group shells are visible within the innermost 1.5 Mq. The bottom panel of Figure [T8] shows 
the mass fraction profiles from the post-process, exponential, and power-law trajectories. Mass 
fraction profiles correspond to the left-hand y-axis, and the peak temperature and peak density 
curves correspond to the right-hand y-axis. Figure|3]also shows the peak temperatures and densities 
of this explosion model, and explains the general trends in the "^'^Ti and ^^Ni mass fraction profiles 



of Figure 18 For interior masses less than « 0.4 Mq, the synthesis of '^'^Ti is due to the or-rich 
and orp-rich freeze-outs, regions 3 and 4 respectively in Figure |4| In this mass range, the "^"^Ti and 
^^Ni mass fractions given by the exponential and power-law profiles have about the same level of 
agreement with the post-process values. Near 0.4 Mq, the thin chasm region separating region 3 
(a-rich freeze-out) and region 6 (silicon-rich) is traversed, which causes the downward spike in 
the "^"^Ti profile (see Figures |3] and |4]). The downward spike occurs at different mass locations for 
the post-process, exponential and power-law profiles because when the thin chasm is traversed 
depends on the exact values of peak temperatures and densities, and the location of the thin chasm 
hne on the temperature-density plane due to the timescale of the expansion profile. 

Between « 0.4 Mq and « 0.7 Mq the explosion is operating in the silicon-rich, region 6 of 
Figure [s] In this region the final yields of "^"^Ti and ^^Ni are sensitive to the initial composition. 
Beyond « 0.7 Mq there are two reasons for the rapid decline of the "^"^Ti and ^^Ni mass fraction 
profiles from the parameterized trajectories while the post-process mass fraction profiles remain 
essentially flat at « 10"^ out to 1.5 Mq First, the parameterized trajectories assume an initial com- 
position that is generally pure ^^Si (with perhaps some neutrons or protons to adjust Y/, see Section 
[3]) which is different than the initial composition of the hypemova model. In region 6 and where 
peak temperature drops below 4x10^ K the yields from the parameterized trajectories are initial 
composition dependent. Second, the strong shock that follows the weak shock raises the tempera- 
ture to ~ 3.5x10^ K at densities of « 10^ g cm~^ out to ~ 1.5 Mq. Over another 5 s of evolution, 
this is sufficient to turn some of the ^^Si into "^"^Ti. 

Yields from the exponential and power-law trajectories are generally only within an order 
of magnitude of the post-process yields. Integrating the '^'^Ti mass fraction profiles in Figure 18 



over the interior mass gives the total mass of "^"^Ti ejected by this model. We find 2.66x10"^ Mq 
for post-processing, 5.23x10^^ Mq for the exponential profile, and 6.34x10"^ Mq for the power- 
law profile. Similarly for ^^Ni, we find 3.77x10"' Mq for post-processing, 4.83x10 ' Mq for 
the exponential profile, and 4.97x10"' Mq for the power-law profile. When the hydrodynamic 
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evolution is complicated, the post-process and parameterized profile yields will generally not agree, 
but parameterized profiles still provide guidance on interpreting the post-processed yields. 



8.3. A 2D Rotating Supernova Model 



Our third model is a 2D explosion of a rotating 1 5 M© star (model 1 from| Fryer & Heger|2000 1 



This simulation follows the collapse, bounce, explosion, and includes the entire proto-neutron star 
throughout the evolution. The dynamical trajectory of the particles, and the time material spends at 
any location, can play a major role in the final "^"^Ti and ^^Ni yields. In some cases, the convective 
motion makes the matter undergo a series of heating-cooling cycles. However, for most matter, 
these cycles occur well above NSE temperatures and, fortunately, the evolution in the NSE regime 
does not affect the final yields significantly except through changes in the electron fraction. This 
rotating supernova model ends at 1.4 s after bounce which has two implications. First, some of 
the nominal ejecta may still fall back onto the proto-neutron star. That is, some of the particles we 
post-process may not ultimately be part of the nucleosynthetic yield. Second, most of the particles 
have temperatures large enough to interfere with a comparison of material that has undergone a 
complete freeze-out via the exponential or power-law trajectory. To facilitate this comparison we 
have appended exponential and power-law thermodynamic trajectories to the final time point of 
the dynamical model. In this manner we extend the thermodynamic evolution to 4.2 s, by which 
time all the particles have temperatures below 5x10^ K. The quantitative differences between the 
exponential tail and power-law tail appear small enough that we will only show results for the 
power-law tail. 



The top row of Figure 19 shows the peak temperatures and peak densities within the innermost 
~ 1 .2 Mq at the coordinates reached by all particles at 4.2 s. The equatorial plane is located at y=0, 
two lobes appear at roughly + 45°, and the overall geometry is not symmetric due to rotation and 
convective fluid motions. Most of the peak thermodynamic conditions are within the bounds of 
our analysis. We ignore those particles with peak temperatures above 1x10'° K or peak densities 
above 1x10^° g cm~^. The post-process yields of '^'^Ti and ^^Ni within the innermost ~ 1.2 Mq are 
shown in the bottom row of Figure 19 Most of the '^'*Ti is created within roughly + 15° of the 



equatorial plane, but split into two distinct regions because of the "^"^Ti chasm, the QSE-leakage 
region 2 of Figure |4j 



The first row of Figure 20 shows the peak temperatures and peak densities as a function of 
interior mass. There is considerable scatter at almost any mass location due to the asymmetries 
inherent in the 2D model. Mass fraction profiles of ''^'^Ti for the post-process, exponential, and 
power-law trajectories are compared in the second row of the figure. Between ^ 0.1 M© and « 0.3 
Mq the "^"^Ti mass fractions from all three thermodynamic trajectories rapidly decrease because a 
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subset of the particles have peak temperatures and peak densities characteristic of the chasm, the 
QSE-leakage region 2 of Figure |4} The width of the chasm associated with an exponential profile 
is narrower than the chasm of a power-law profile, accounting for the power-law "^"^Ti mass fraction 
profile decreasing more rapidly than the exponential "^"^Ti mass fraction profile in this mass range. 
As expected from our previous analysis, ^^Ni undergoes no such phase transition, with the result 
that all thermodynamic trajectories give values of order unity in this mass range. Between « 0.6 
Mq and « 0.7 Mq a large subset of the particles in the equatorial regions have peak temperatures 



that drop below 4x10 K (see Figure 20) and traverse the thin chasm. This results in the strong 



trend towards decreasing "^"^Ti and ^^Ni mass fraction profiles in this mass range. Beyond « 0.6 



Mq the two asymmetric lobes of Figure 19 are visible as the two horizonal bands in the second 
row of Figure |20j The "rain" of points in this mass range is due to the large scatter in the peak 
thermodynamic conditions, where a number of particles have peak temperatures that drop below 
4x10^ K and traverse the thin chasm. 

Integrating the '^'^Ti mass fraction profiles over the interior mass gives the total mass of "^"^Ti 
ejected by this model. We find 6.98x10^^ Mq for post-processing, 5.09x10"^ Mq for the exponen- 
tial profile, and 4.82x10"^ Mq for the power-law profile. Similarly for ^""Ni, we find 3.89x10"^ Mq 
for post-processing, 3.99x10"' Mq for the exponential profile, and 4.10x10"' Mq for the power- 
law profile. 

Overall, the mass fraction profiles from exponential and power-law trajectories are generally 
within a factor ~ 4 of the post-process values, except in regions where the chasm is traversed, and 
the total yields of ''''Ti and ^^Ni from the parameterized profiles mimic the post-process yields for 
this 2D, rotating supernova model. 



9. Discussion 

We have explored the trends in, and sensitivity to, the '^'^Ti and ^^Ni yields in the ejecta of three 
contrasting core-collapse supemovae models. We used yields from two parameterized expansion 
profiles and compared them to the yields from post-processing trajectories from the core-collapse 
models. 

Both parameterized profiles, the classic exponential and our new power-law expressions, as- 
sume a constant Ip adiabat throughout the evolution. For any given peak temperature and peak 
density initial conditions, the power-law is slower compared to the exponential and together they 
generally bound the trajectories from core-collapse simulations. We find that ''"'Ti may be produced 
by more than one type of freeze-out, depending on the peak temperatures, densities and electron 
fraction values of the thermodynamic trajectories. We have identified several distinct regions in 
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the peak temperature-density plane from the parameterized profiles. Each region is characterized 
by different types of transitions that the QSE cluster sustains during the evolution. Reactions 
that break equilibrium are responsible for the flow transfer to the remaining small scale clusters, 
maintaining their structure until freeze-out. The result is unique mass fraction profiles per region. 
The transitions are entropy driven, not expansion timescale driven, where the temperature sets an 
approximate threshold for a transition, while the density at the threshold temperature determines 
whether the transition takes place or not. The expansion timescale affects the locus of the borders 
among different regions in the peak temperature-density plane. 

For initially symmetric matter, region 1 is the normal freeze-out regime, where no phase tran- 
sition takes place and the yields from the parameterized profiles are in good agreement with NSE 
or QSE estimates. Region 2 is the ''^Ti chasm, where '^'^Ti is depleted as a result of the large scale 
QSE cluster dissolution to two smaller ones, and the subsequent flow leakage from the Si-group 
towards the Fe-group nuclei. Region 3 is the or-rich freeze-out regime, where the large scale QSE 
cluster shrinks and shifts upward in mass due to the domination of the 3a forward rate over its in- 
verse. Region 4 is the ap-rich freeze-out regime, where p(e~, Ve)n and n(e+, v^)p drive the material 
slightly proton-rich early in the evolution. Region 5 is the regime where photodisintegrations dom- 
inate capture reactions. Region 6 is the incomplete Si-burning regime, where the timescales for the 
plasma to reach large scale QSE or NSE are comparable to the freeze-out timescale, preventing in 
general their establishment. 

Compared to the symmetric case, "^"^Ti and ^^Ni are gradually underproduced for initial < 
0.5 and become less sensitive to reaction rates. The basic structure of the temperature-density plane 
is maintained, with the exception that the ap-rich freeze-out region decreases in size until complete 
extinction for Yg » 0.46 due to the absence of free protons in neutron-rich environments, and the 
"^Ti chasm region expands as a result of the ''^Ti underproduction in the neighboring regions in the 
temperature-density plane. In addition, the region of an-rich freeze-out appears, where both '^'^Ti 
and ^^Ni are depleted since they are symmetric isotopes. 

For initial Y^ > 0.5 '^'^Ti and ^^Ni are favored by large scale NSE and QSE equilibria due to the 
minimization of the Helmholtz energy. This results in ^^Ni still being one of the dominant yields, 
although symmetric isotopes are not favored by the major flows once the large scale QSE cluster 
dissolves. In this electron fraction regime, weak interactions are crucial to '^'^Ti production. The 
dominance of n(e"^, v^)p over p(e , Ve)n early in the evolution and the rest of the weak interactions 
later on, results in significant production of symmetric isotopes. Regions 3 and 4 merge to become 
a regime of ap-rich freeze-out due to the large proton excess, and the '^'^Ti chasm expands. Region 
1 is no longer a normal freeze-out regime, but it is characterized by a phase transition due to the 
interplay between (p, y) and weak reactions. 

The three core-collapse models we post-processed were a ID Cassioppeia A model, a ID 
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double-shock hypernova model, and a 2D rotating 15 Mq model. Mass fractions of '^'^Ti and ^^Ni 
from the exponential and power-law trajectories were shown to generally lie within a factor ~ 
8 or less of the post-process yields, except in regions where the thin chasm is being crossed or 
the temperature fell below 4x10^ K. The total ejected masses of '^'^Ti and ^^Ni were shown to be 
within a factor of ~ 2 or less for all three models. When the thermodynamic trajectories of a core- 
collapse model have an expansion profile similar to a parameterized expansion of any form, it is 
generally safe to trust the yields from the parameterized profiles. For more complicated thermo- 
dynamic trajectories, the yields from the parameterized profiles should not be trusted, although the 
parameterized profiles may provide useful information about the underlying physics. 

The location of the "^"^Ti chasm region is profile dependent, and its width is minimized for 
initially symmetric matter but expand dramatically for These trends could account in 

part for the observed paucity of supernova detected in the light of radioactive "^"^Ti. A mass cut in 
supemovae models where the electron fraction begins to fall below 0.5 may not be the most suitable 
choice, since the layers above the mass cut are biased to initially symmetric compositions and can 
produce ample "^"^Ti. We find that variations in the "^^^Ti/^^Ni ratio originate from variations in '^'^Ti, 
since ^^Ni is produced in large quantities over most of the peak temperature-density plane (Figure 



1 1 1. Although the massive production of ^^Ni and its decay to ^^Co and ^^Fe outshine every other 
decay during their lifetime, the decay of '^'^Ti to '^'^Sc and '^'^Ca has a longer lifetime. This implies 
that measurements of the yield from "^"^Ti may be used to estimate the yield of ^^Ni (e.g. measure 
"^"^Ti in Cas A and deduce Cas A's ^^Ni yield), assuming the supernova models' thermodynamic 
trajectories approximate one of the parameterized profiles and '*'^Ti originates from regions where 
^^Ni was dominant. 



Woosley et al.| ( |1973| ) first defined the notions of the normal and or-rich freeze-outs. Using 



an exponential profile they describe how the various types of freeze-out are driven by the effect 
of the 3a reaction. They identify 3 regions in the peak temperature-density plane - normal and 
a-rich freeze-outs and the incomplete silicon burning regime - based on the availability of a- 
particles. These 3 regions are bordered with thin lines based on semi-analytical relationships for 
the timescale required to reach QSE. Their description is constrained to the interplay between the 
3a reaction and the large-scale QSE cluster. 



Woosley & Hoffman| ( [T992[ ) explored cases of freeze-outs starting at large neutron excesses, 
resulting in the production of isotopes heavier than the Fe-group and illustrating the smooth merg- 
ing of the freeze-out process into the r-process. They used the exponential profile for a peak tem- 
perature = 10, three different values for the peak density and a grid of initial neutron excesses 
within the range < ?/ < 0.21 (corresponding to a range for the electron fraction 0.395 < < 0.5). 
For the specific case of initial Ye = 0.48 they mentioned that yields for the Fe-group nuclei were 
not very different from the symmetric case, but the production of heavier nuclei is also possible. 
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This is in agreement with our results (Table [2]), where the yields for '^'^Ti and ^^Ni are decreased 
but are not negligible compared to the symmetric case, and a region of om-rich freeze-out appears 
in the peak temperature-density planes. Their argument for flows being inhibited beyond ^^Ni for 
nearly symmetric compositions is functionally the same as ours. Their argument is based on Q val- 
ues and other nuclear structure effects around the doubly magic nucleus ^^Ni, which is equivalent 
to a persistent small equilibrium cluster localized around ^^Ni until complete freeze-out. Finally, 
they suggested the possibility of freeze-out expansions with negative neutron excesses (7^ > 0.5) 
due to the interplay between p(e~,Ve)n and n(e^,v^)p, which we have discussed in detail within 
this study. 



The et al. ( 1998 1 conducted the first detailed sensitivity study of reaction rates on the "^"^Ti 
mass fraction by using the exponential profile and one pair of peak conditions (Tg = 5.5, p = 10^ 
g cm~^). They identified a significant amount of reactions affecting "^"^Ti, and suggested the impact 
of "^^V(p, y)'*^Cr and ^^C{a,yy^O for Ye < 0.5. However, the choice of only one set pair of peak 
conditions and one expansion profile doesn't allow all crucial reactions identified in all regions. 



Recent sensitivity studies have varied rates within their experimental uncertainty limits. Hof- 



fman et al. (2010) used one-zone calculations with results for '^'^Ti, ^^Ni, ^^Ni and ^^Ni presented 
in proportion to solar ^^Fe, considering in addition yields for these isotopes from published super- 
novae models. They varied all published rate compilations of '*'^Ti(ar, p)'^^V and ^^Ca{a, y)'^"^Ti for 
peak conditions taken from various points within the temperature-density plane given in (Magkot^ 
sios et al.] ( 2008[ ). They also varied the expansion timescale of the exponential profile, concluding 
that the yield of "^"^Ti is affected by the time it remains within successive burning stages. Tur et al. 



(20101) have utilized stellar evolution models to address the impact of the 3a and ^^C(Qf, y)^^0 
reactions on ^^Al, '^'^Ti and ^"Fe. They conclude that "^"^Ti is relatively insensitive to these rates, 
although they comment that this result depends on the explosion physics and supernovae rate as- 
sumed. Variations to the 3a rate within our study also result in relatively insensitive yields, because 
only a limited flow is required by this reaction to impact the QSE cluster. In order to illustrate the 
reaction's role, it is necessary to remove it completely from network calculations. ^^C{a, yY^O has 
an impact on '*'*Ti mostly for neutron-rich compositions (Table [s]). 

Our study adds to these works by considering all the freeze-out regions in the peak temperature- 
density plane over a broad range of Ye. The various types of freeze-out may be understood as dif- 
fering equilibrium patterns during the evolution, where a change to the pattern is usually signaled 
by individual reactions dropping out of equilibrium. The use of two expansion profiles and detailed 
sensitivity studies for all regions within the peak temperature-density plane reveals the importance 



of additional reactions crucial to "^"^Ti synthesis beyond the set first identified by The et al. ( 1998 1. 
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Tablel. NUCLEAR NETWORKS 
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A ■ A A ■ A A ■ A A ■ A 
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Ga 


61 


70 


59 


75 


59 


94 


56 


92 


Ge 


64 


71 


62 


78 


62 


97 


58 


95 


As 






65 


79 


68 


103 


60 


100 


Se 






67 


83 


63 


103 


63 


103 


Br 






68 


83 


69 


106 


65 


105 


Kr 






69 


87 


72 


109 


68 


108 


Rb 






73 


85 


74 


113 


74 


111 


Sr 






74 


84 


76 


118 


73 


114 


Y 






75 


87 


78 


121 


75 


119 


Zr 






78 


90 


80 


122 


77 


122 


Nb 






82 


90 


81 


123 


80 


124 


Mo 






83 


90 


82 


125 


82 


125 


Tc 






89 


91 


87 


127 


85 


126 


Ru 










90 


130 


86 


128 


Rh 










93 


131 


89 


130 


Pd 










94 


132 


91 


132 


Ag 










97 


133 


93 


134 


Cd 










98 


136 


95 


134 


In 










99 


149 


97 


137 


Sn 










102 


152 


99 


140 


Sb 














105 


151 


Te 














109 


154 


I 














112 


158 


Xe 














113 


161 


Cs 














118 


165 


Ba 














119 


168 


La 














122 


172 


Ce 














122 


175 


Pr 














126 


178 


Nd 














127 


178 



Table 1 — Continued 



204 489 1341 3304 





^max ^min ^max ^min ^max ^min 


^max 


Pm 


no 


185 


Sm 


133 


188 


Fu 


136 


190 


Gd 


139 


191 


Tb 


142 


192 

J- ^ ^ 


Dv 


1 43 


193 


Ho 


146 


196 

J- y\J 


Er 


148 


198 


Tm 


150 


198 


Yb 


152 


200 


T u 


156 


209 


Hf ••• 


159 


212 


Ta ••• 


161 


217 


W 


163 


220 


Re ••• 


167 


225 


Os ••• 


169 


226 


Ir 


172 


230 


Pt 


175 


232 


Au ••• 


178 


236 


Hg ••• 


178 


239 


Tl 


182 


245 


Pb ••• 


185 


246 


Bi 


188 


251 


Po ••• 


193 


237 


At ••• 


210 


211 
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Table 2. FREEZE-OUT DOMINANT YIELDS 



Freeze-out 


Ppeafc (g cm ^) 


Exponential 


Power-Law 


Ye = 0.48, = 9 


Normal 


5x 10'^ 


^^Fe, ^^Ni, ^^Fe, ^^Fe, *^°Ni 


^^Fe, 5^Fe, ^^Ni, '^"Ni, ^^Fe 


a-rich 




5^Ni, ^He, ''"Ni, ^^Zn, ^^Zn 


5«Ni, ^He, ''"Ni, '"^Zn, ^^Zn 


ap-rich 








orn-rich 


5 X 10^ (79 = 6) 


^He, «'Kr, n, ^^Ca, ^^Se 


^He, n, 49Ca, ^^Kr, '^Zn 


Y, = 0.50, = 9 


Normal 


5x 10'' 


5^Ni, ^^Fe, ^^Fe, ^^Co, "Ni 


54Fe, ^**Ni, 5^Ni, ^-^Fe, ^^Fe 


a-rich 


10« 


^^Ni, ^He, 6°Zn, "Ni, ^^Cu 


5^Ni, ^^Cu, ^He, "Ni, ^^Ni 


ap-rich 


10^ 


^He, ^^Ni, p, ^°Zn, ^^Pe 


^He, 5^Ni, p, 60Cu, ^^Fe 


an-vich 








Ye = 0.52, = 9 


(p, 7) leakage 


10^ (Tg = 5) 


56Ni, 5^Cu, ^'^Cu, 58Ni, ^^Ni 


5^Ni, ^'^Cu, ^''Ni, ^°Zn, ^^Ni 


or-rich 








orp-rich 


10' 


^^Ni, ^He, p, ^°Zn, ^^Ni 


^^Ni, ^He, p, 5'Ni, ^Cu 



orn-rich 
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Table 3. NUCLEAR REACTIONS RELEVANT TO ^Ti SYNTHESIS 



Reaction 


Contribution 


Rank 


Ye 


Region 


Profile 


Global scope reactions 


3a 


flow transfer to QSE cluster 


primary 


0.48-0.52 


2-5 


both 


p(e",Ve)n/n(e+,v;)p 


Ye adjustment 


primary 


0.48-0.52 


1-5 


both 


p(e",Ve)n/n(e+,y^)p 


chasm widening 


primary 


0.50-0.52 


2 


both 


(a, y) reactions 


^''Ca(Qr,y)^^Ti 


2nd arc amplitude/slope 


primary 


0.48-0.52 


3-4 


both 


i2C(Qr, 7)1^0 


flow transfer to QSE cluster 


secondary 


0.48 


2-5 


both 


'Be(a,r)>iC 


chasm depth 


secondary 


0.5 


2 


both 


24Mg(a,y)28Si 


chasm depth 


secondary 


0.5 


2 


both 


^2Ca(a,y)4^Ti 


chasm depth, 1st arc dip 


secondary 


0.48 


2-3 


both 


(a, p) reactions 


^4Ti(Q',p)'*^V 


chasm formation, depth, shift 


primary 


0.48-0.52 


2 


both 




1st arc dip / 2nd arc slope 


primary 


0.48-0.50 


3-4 


both 


4^Ti(a,p)4^V 


1st arc dip 


primary 


0.52 


1 


both 




2nd arc dip/slope 


primary 


0.48-0.50 


3-4 


both 


i^F(Qr,p)20Ne 


1st arc dip 


primary 


0.52 


1 


both 


2iNa(a,p)24Mg 


1st arc dip 


primary 


0.52 


1 


both 


^'^Ca(Qr,p)43sc 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


2^Al(a,p)3"Si 


chasm depth 


secondary 


0.5 


2 


exponential 


55Co(Q',p)^**Ni 


chasm depth 


secondary 


0.5 


2 


exponential 


^^Cr(a,p)^iMn 


chasm depth 


secondary 


0.5 


2 


exponential 


i^F(af,p)20Ne 


chasm depth 


secondary 


0.5 


2 


exponential 




chasm depth 


secondary 


0.5 


2 


exponential 


54Fe(Qr,p)57Co 


chasm depth 


secondary 


0.5 


2 


exponential 


2iNa(Qr,p)24Mg 


chasm depth 


secondary 


0.5 


2 


exponential 


56Ni(ar,p)59Cu 


chasm depth 


secondary 


0.5 


2 


exponential 


^Li(a,p)''Be 


chasm depth 


secondary 


0.5 


2 


power-law 


"N(a,p)i60 


chasm depth 


secondary 


0.5 


2 


power-law 


^2Ca(a,p)45sc 


1st arc dip / 2nd arc slope 


secondary 


0.48 


3 


both 


^3Sc(a,p)46Ti 


1st arc dip 


secondary 


0.48 


3 


both 


^^Ni(Qr,p)^iCu 


2nd arc amplitude 


secondary 


0.48 


3 


both 



-45- 



Table 3 — Continued 



Reaction 


Contribution 


Rank 




Region 


Profile 


38Ca(a,p)^iSc 


1st arc dip 


secondary 


0.52 


1 


both 


34Ar(Qr,p)37K 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


3^Ca((^,p)^iSc 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


(p, 7) reactions 


^^V(p,7)^^Cr 


2nd arc formation/dip 


primary 


0.50-0.52 


1,3,4 


both 




3rd arc formation 


primary 


0.52 


4 


power-law 


^'Sc(p,y)42Ti 


2nd arc dip 


primary 


0.5 


4 


both 


^^Sc(p, yf^Tx 


2nd arc dip 


primary 


0.5 


4 


both 


^^Ti(p,y)45v 


2nd arc dip / 3rd arc formation 


primary 


0.5 


4 


both 


'*^Ti(p,y)45V 


44xi _ cluster 


primary 


0.5 


4 


both 


4iSc(p,r)42Ti 


2nd arc dip 


primary 


0.5, 0.52 


4 


both 


^^Ni(p,7)5^Cu 


flow transfer within QSE cluster 


primary 


0.5 


3 


both 


^5V(p,7)4^Cr 


regions 3-4 borderline 


secondary 


0.5 


3-4 


both 


^'^Ca(p,7)4iSc 


post-2nd dip track 


secondary 


0.5 


4 


both 


'*^V(p,7)45Cr 


post- 2nd dip track 


secondary 


0.52 


1,4 


both 


43Ti(p,7)44v 


post-2nd dip track 


secondary 


0.52 


1,4 


both 


42Sc(p,7)«Ti 


post- 2nd dip track 


secondary 


0.52 


1,4 


both 


5^Cu(p, 7)^^Zn 


post-2nd dip track 


secondary 


0.52 


1,4 


both 


20Ne(p,7)2iNa 


post-2nd dip track 


secondary 


0.52 


1,4 


both 


^'V(p,7)4^Cr 


post-2nd dip track 


secondary 


0.52 


1,4 


both 


43Sc(p,7rTi 


post-2nd dip track 


secondary 


0.48, 0.52 


1,3,4 


both 


«Ti(p,7rv 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


4iSc(p, 7)42^1 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


^3Sc(p,7)44Ti 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


^«Ca(p,7)4iSc 


post-2nd dip track 


secondary 


0.52 


4 


exponential 


^«Ca(p,7)4iSc 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


^2Ca(p,7)43sc 


2nd arc slope 


secondary 


0.48 


3 


both 


39K(p,7)4"Ca 


2nd arc slope 


secondary 


0.48 


3 


exponential 


^^Co(p,7)5^Ni 


chasm depth 


secondary 


0.5 


2 


exponential 


54Fe(p,7)55Co 


chasm depth 


secondary 


0.5 


2 


exponential 


52Mn(p,7)53Fe 


chasm depth 


secondary 


0.5 


2 


exponential 
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Table 3 — Continued 



Reaction 


Contribution 


Rank 




Region 


Profile 


5^Ni(p,y)5^Cu 


regions 3-4 borderline 


secondary 


0.5 


3-4 


both 


weak reactions 




,Ve)5^Fe 


chasm widening 


primary 


0.50-0.52 


2 


both 


5°Mn(e" 


",Ve)^°Cr 


chasm widening 


primary 


0.50-0.52 


2 


both 


^^Ni(e- 


,Ve)^^Co 


chasm widening 


primary 


0.50-0.52 


2 


both 


5^Cu(e- 


,Ve)^'Ni 


chasm widening 


primary 


0.50-0.52 


2 


both 


53Co(e- 


, Ve)"Fe 


chasm widening 


primary 


0.50-0.52 


2 


both 


5^Fe(e- 


Ve)''Mn 


chasm widening 


primary 


0.50-0.52 


2 


both 




,Ve)5^Ni 


chasm widening 


primary 


0.50-0.52 


2 


power-law 


^^Cu(e- 


,Ve)5^Ni 


chasm widening 


primary 


0.50-0.52 


2 


power-law 


^2Ti(e- 


,Ve)^'Sc 


2nd arc dip / post-2nd dip track 


primary 


0.50-0.52 


4 


both 


^iSc(e- 


,Ve)^'Ca 


2nd arc dip / post-2nd dip track 


primary 


0.50-0.52 


4 


both 


^3Ti(e- 


,Ve)^^Sc 


2nd arc dip / post-2nd dip track 


primary 


0.50-0.52 


4 


both 




,Ve)'^Ti 


2nd arc dip / post-2nd dip track 


primary 


0.50-0.52 


4 


both 




,Ve)'"Ti 


2nd arc dip / post- 2nd dip track 


primary 


0.52 


4 


exponential 


3^Ca(e" 




2nd arc dip / post- 2nd dip track 


primary 


0.52 


4 


power-law 


3''Ca(e" 


",Ve)''K 


2nd arc dip / post- 2nd dip track 


primary 


0.52 


4 


power-law 


3"S(e- 


,Ve)^"P 


flow transfer to symmetric nuclei 


primary 


0.52 




both 


5^Cu(e- 


,Ve)5^Ni 


flow transfer to symmetric nuclei 


primary 


0.52 




both 


34Ar(e- 


,ye)''Cl 


flow transfer to symmetric nuclei 


primary 


0.52 




both 


42Ti(e- 


,ye)''Sc 


flow transfer to symmetric nuclei 


primary 


0.52 




both 


4^V(e- 


,Ve)^^Ti 


flow transfer to symmetric nuclei 


primary 


0.52 




exponential 


^5Cr(e- 




flow transfer to symmetric nuclei 


primary 


0.52 




power-law 


55Ni(e- 


,Ve)5^C0 


flow transfer to symmetric nuclei 


secondary 


0.52 




power-law 


58Zn(e- 


,Ve)^'Cu 


flow transfer to symmetric nuclei 


secondary 


0.52 




both 


43Ti(e- 




flow transfer to symmetric nuclei 


secondary 


0.52 




power-law 


4iTi(e- 


,Ve)^lSc 


flow transfer to symmetric nuclei 


secondary 


0.52 




power-law 


45v(e- 


,Ve)^^Ti 


flow transfer to symmetric nuclei 


secondary 


0.52 




power-law 


3^Ca(e- 




flow transfer to symmetric nuclei 


secondary 


0.52 




power-law 


3^Ca(e- 




2nd arc dip / post-2nd dip track 


secondary 


0.52 


4 


power-law 


34Ar(e- 


,ye)''Cl 


2nd arc dip / post-2nd dip track 


secondary 


0.52 


4 


power-law 
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Table 3 — Continued 



Reaction 


Contribution 


Rank 




Region 


Profile 


5'Cu(e~,ye)"Ni 


2nd arc dip / post- 2nd dip track 


secondary 


0.52 


4 


power-law 


3^Ar(e-,ye)3^CI 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


3^Ca(e-,Ve)^^K 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


3^Ca(e-,Ve)^^K 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


^2Ti(e-,Ve)^2Sc 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


5^Cu(e-,Ve)^^Ni 


3rd arc amplitude 


secondary 


0.52 


4 


power-law 


(p, n) and (a, n) reactions 


5^Co(p, n)5^Ni 


chasm widening 


secondary 


0.5 


2 


both 


5^Co(p, n)56Ni 


chasm widening 


secondary 


0.5 


2 


both 


2^Al(p,n)27Si 


chasm widening 


secondary 


0.5 


2 


both 


"B(p,n)iiC 


chasm widening 


secondary 


0.5 


2 


both 


iOB(Qr,n)i3N 


chasm widening 


secondary 


0.5 


2 


both 


"B(a,n)i4N 


chasm widening 


secondary 


0.5 


2 


both 


20Ne(Qr,n)23Mg 


chasm widening 


secondary 


0.5 


2 


both 


9Be(a, ny^C 


chasm widening 


secondary 


0.5 


2 


both 


^2Ca(a,n)^5Ti 


1st arc dip / 2nd arc slope 


secondary 


0.48 


3 


both 


34S(a,n)3'Ar 


2nd arc ampUtude 


secondary 


0.48 


3 


both 
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Fig. 1. — Reaction Q values for the (or, y) channels within our reference network of 489 isotopes. 
Qay is also equivalent to the alpha capture thresholds Sa- (N,Z) boxes correspond to the Q value 
of reaction "^Xia, j)^^ ■ White boxes imply the absence of an (or, y) reaction. 
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Fig. 2. — Schematic temperature or density evolution for the exponential (dashed) and power-law 
(solid) profiles. For either of the two initial values illustrated, the exponential profile declines faster 
than the power-law profile. Passages through different burning regimes are indicated. 
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Peok Temperoture (T, 



Peok Temperoture (Tg) 



Fig. 3. — Final mass fraction of '^'^Ti (first row), ^^Ni (second row) and '^He (third row) in the peak 
temperature-density plane for the exponential thermodynamic profile (first column) and power- 
law profile (second column) at Ye=0.5. Different colored triangles show the temperature-density 



positions from different supernova and hypernova models - blue for a ID Cas A model (Young 
et al.|2006 ), gray for the 2D rotating progenitor E15B model ( Fryer & Heger|2000 ), pink for a ID 
hypernova model ( Fryer et al.|2006 ), and cyan for a 2D magnetohydrodynamic coUapsar model. 
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Peok Temperoture (Tg) 



Fig. 4. — Final mass fraction of Ti in the peak temperature-density plane for the exponential 
thermodynamic profile at Ye=0.5. Six distinct regions of "^"^Ti synthesis are labeled. Region 1: 
normal freeze-out from NSE, abundance largely determined from Q values. Region 2: Chasm 
region, passage from 1 QSE cluster to 2 QSE clusters. Region 3: a-rich freeze-out. Region 4: ap- 
rich freeze-out. Region 5: Photodisintegration regime, neutrons, protons, and a dominate. Region 
6: Incomplete silicon burning, ^^Si rich. The thin cyan line is the locus of points where tqse = 
0.012 

^freeze • 




^Sw; — . ^hi 




„ 10" 



8 7 6 


5 4 3 2 1 C 

Tg 






/.i4 vy^ 

t \ 



876543210 



7 6 5 4 3 2 1 



Fig. 5. — Mass fraction evolution of ^*^Si, "^"^Ti, ^^Ni, neutrons, protons, and a-particles for the 
exponential (left column) and power-law (right column) expansions at ¥^=0.5. In each panel the 
peak temperature is 9x10^ K. The top row corresponds to region 1 in Figure |4](p = 5 x 10^ g 
cm"^), second row to region 2 (p = 1.75 x 10^ g cm"^), third row to region 3 (p = 1 x 10^ g cm"^), 
and fourth row to region 4(p = 1x1 O^g cm"^). 
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Fig. 6. — Normalized nuclear flows at three different points in the evolution of "^^^Ti from freeze- 
out; before the minimum (top panel), at the minimum (middle panel), and after the minimum 
(bottom panel). Normalized flows ^ are colored black for < ^ < 0.01, navy for 0.01 < < 0.05, 
blue for 0.05 < < 0.1, cyan for 0.1 < (p < 0.4, green for 0.4 < < 0.8, red for 0.8 < < 1.0, 
yellow for 0=1.0. Normalized flows in the vertical direction, corresponding to {p, y) reactions, and 
along diagonal directions corresponding to (a, p) reactions have been drawn thicker for clarity. 




Fig. 7. — Reaction network (solid curves) and QSE solutions (dashed and dot-dashed curves) for 
'^'^Ti (upper left), "^^Ca (lower left), and "^^Cr (lower right). Each element shows the characteristic 
arc driven by equilibrium (a, p) and (p, y) reactions. Beginning at the local minima and onwards 
to lower temperatures only (p, y) reactions remain in equilibrium, which drive the mass fractions 
to larger values. This general behavior applies to most of the elements within the silicon and iron 
groups (small sample upper right). 



-55- 




Fig. 8. — Mass fraction evolution of key isotopes from network and NSE calculations for peak 
temperatures and peak densities corresponding to a normal freeze-out (top panel) and an a-rich 
freeze-out (bottom panel). 
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Fig. 9. — Final yield of Ti in the peak temperature-density plane for different values of the 
initial electron fraction Ye under the exponential freeze-out profile. The top row, from left to right, 
corresponds to Ye = 0.484, 0.486, and 0.488. The second row, from left to right, corresponds to Ye 
= 0.490, 0.492, and 0.494. The third row, from left to right, corresponds to Ye = 0.496, 0.498, and 
0.500. The bottom row, from left to right, corresponds to Ye = 0.502, 0.504, and 0.506. 
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Fig. 10. — Same as Figure |9} but for the final yield of ^^Ni. 




Fig. 11. — Same as Figure |9} but for the "^"^Ti normalized production factor P44. 




Fig. 12. — Final mass fraction of ^"^Ti (first row), ^^Ni (second row) and "^He (third row) in the peak 
temperature-density plane for the exponential thermodynamic profile at yg=0.48 (left column) and 
7^=0.52 (right column). The distinct regions of "^"^Ti synthesis are labeled. Region 1: normal 
freeze-out from NSE (left), (p,y) leakage from symmetric to proton-rich nuclei (right). Region 2: 
Chasm region, passage from 1 QSE cluster to 2 QSE clusters. Region 3: a-rich freeze-out. Region 
4: ap-rich freeze-out. Region 5: Photodisintegration regime, neutrons, protons, and a dominate. 
Region 6: Incomplete silicon burning, ^^Si rich. Region 7: an-rich freeze-out. 
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Fig. 13. — Mass fraction of select isotopes in NSE at = 9 and p = 10^ g cm ^ for different 
electron fractions. 
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Fig. 14. — Examples of Ti sensitivity to reaction rates for the power-law thermodynamic profile 
at 7^=0.5. From left to right in the first row, the contour plots show the effects of removing the 
3a and increasing the p(e~,ye)n + n(e'^,v^)p by 1000, while the third plot shows the ratio of the 
'^'^Ti yield with single rates increased by 100 to the nominal '^'^Ti yield during or-rich freeze-out. 
Isotopes whose reaction rates produce interesting variations are labeled. Black color is used for 
the (a,y) reactions, green for the {a,p) reactions, red for (p,y) reactions, cyan for (a,n), and blue 
for the weak reactions. The second row shows from left to right the '^'^Ti mass fraction sensitivity 
to '*'*Ti(Q:, p)'^^V and "^"CaCa, y)'*'^Ti for the or-rich freeze-out, and '^'^Ca(a,y)'^'^Ti for the ap-hch 
freeze-out. Black curves are for the nominal rates, red curves are for the QSE yields, green curves 
for rates multiplied by 0.01 and blue curves for rates multiplied by 100. 




Fig. 15. — Evolution of the electron fraction versus temperature (upper left) and time (up- 
per right) for the exponential and power-law profiles starting from conditions, Tg = 10 GK, 
p = 10'° g cm"^ and = 0.5, corresponding to the normal freeze-out regime (region 1). The 
corresponding evolution of the mass fractions are shown for the exponential profile (lower left) 
and power-law profile (lower right). 
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Fig. 16. — Evolution of the electron fraction for the exponential and power-law freeze-out 
profiles for a freeze-out starting from Tg = 10 GK, p = 10'" g cm"^ and = 0.52, corresponding 
to the (p,7) leakage regime (region 1). The evolution is shown versus temperature (upper plot) and 
time (lower plot). 
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Fig. 17. — Mass fractions from post-processing the Lagrangian thermodynamic trajectories (top 
panel) and post-processing vs parameterized profiles (bottom panel) for the 16 Mq Cas A model 
( Young et al.|2008 1. See text for a discussion of the sharp dip in "^"^Ti « 0.2 M©, and the decline of 
^^Ti and ^^Ni « 0.6 Mq. 
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Fig. 18. — Mass fractions from post-processing the Lagrangian thermodynamic trajectories (top 
panel) and post-processing vs parameterized profiles (bottom panel) for the spherically symmetric 
40 Mq hypernova model ( jFryer et al.|2006[ ) that features a weak shock followed by a strong shock, 
with our reference 489 isotope network. Note the sharp dip in "^"^Ti « 0.4 Mq, and the extended '^'^Ti 
and ^^Ni distributions for the post-processed trajectories. 
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Fig. 19. — Peak temperatures and peak densities (top row) and mass fractions of "^"^Ti and ^^Ni 



(bottom row) generated by the innermost regions in a 2D explosion of a rotating 15 Mq star (Fryer 



& Heger 2000). Particles coordinates are shown at 1.4 s, the end of the dynamical model. Note 



the y=0 equatorial plane, the double lobbed structure, the general asymmetry of the model and 
evidence for the '^'^Ti chasm along the equatorial regions. 
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Fig. 20. — Peak temperatures in grey and peak densities in purple (first row), and comparison 
of post-process vs parameterized "^"^Ti (second row) and ^^Ni (third row) profiles as a function of 
interior mass for the 2D explosion of a rotating 15 Mq star. The scatter is due to the asymmetries 
of the model. 



